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This  dissertation  is  concerned  with  the  development  of  new  methodology  for 
dynamic  lot  size  models. 

We  propose  a new  forward  dynamic  programming  algorithm  with  continuous 
state  variable  (cumulative  production)  for  dynamic  lot  size  models  (DLSM)  in  which 
production  and  inventory  cost  functions  are  only  assumed  to  be  piecewise  linear.  In 
particular,  there  are  no  assumptions  of  convexity,  concavity  or  monotonicity.  Ar- 
bitrary capacities  on  both  production  and  inventory  may  occur,  and  backlogging  is 
allowed.  Thus  the  algorithm  addresses  most  variants  of  the  DLSM  appearing  in  the 
literature.  Computational  experience  shows  it  to  be  very  effective  on  NP-hard  ver- 
sions of  the  problem.  For  example,  48-period  capacitated  problems  with  production 
costs  defined  by  eight  linear  segments  are  solvable  in  less  than  2.5  minutes  of  Vax 
8600  cpu  time. 

We  then  extend  the  algorithm  for  models  with  two  state  variables.  The  new 
state  variable  is  related  to  machine  status,  for  example,  machine  age  or  total  number 
of  setups  or  whether  the  machine  is  on  or  not.  Thus,  the  DLSM  with  learning  in 
setups,  the  DLSM  with  startup  and  reservation  costs,  as  well  as  the  DLSM  with 
machine  replacement  can  be  solved  more  effectively  than  with  existing  algorithms  in 
the  literature. 
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This  dissertation  also  investigates  the  issue  of  the  error  bound  of  the  DLSM. 
Unlike  the  literature  that  computes  the  error  bound  after  the  production  plan  of 
the  first  few  periods  is  determined,  our  approach  decides  the  plan  for  which  the 
corresponding  error  bound  is  the  smallest. 

As  a natural  extension,  we  propose  a heuristic  algorithm,  based  on  the  algorithm 
for  the  single-item  DLSM,  for  the  multi-item  DLSM.  Preliminary  computational  ex- 
perience shows  that  this  algorithm  is  effective  compared  to  existing  algorithms. 
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CHAPTER  1 
INTRODUCTION 


This  research  is  concerned  with  the  development  of  new  methodology  for  dy- 
namic lot  size  models. 

The  contribution  of  this  research  is  summarized  by  the  following: 

1.  A new  dynamic  programming  algorithm  is  developed  for  the  single-item  capac- 
itated dynamic  lot  size  model.  This  algorithm  can  solve  the  single-item  lot  size 
problems  with  hundreds  of  periods  in  reasonable  computational  effort.  In  the 
existing  literature,  these  problems  are  intractable  when  the  number  of  periods 
is  more  than  one  hundred. 

2.  This  research  also  extends  the  idea  of  the  algorithm  for  the  lot  size  model  with 
two  state  variables,  cumulative  production  and  machine  status.  As  a result,  the 
lot  size  models  with  learning  in  setups  and  models  with  startup  and  reservation 

costs,  as  well  as  models  with  machine  replacement,  can  also  be  solved  very 
efficiently. 

3.  For  the  error  bound  of  the  single-item  uncapacitated  dynamic  lot  size  model,  this 
research  provides  a new  algorithm  that  can  obtain  the  error  bound  in  0(T  log  T ) 
where  T is  the  number  of  periods.  In  the  literature,  existing  algorithms  can 
compute  only  the  error  bound  for  the  problems  where  the  unit  production  in  a 
period  is  always  not  greater  than  the  summation  of  the  unit  production  and  the 
unit  holding  cost  in  the  previous  period.  Moreover,  this  research  investigates 
how  a production  plan  up  to  the  decision  horizon  is  chosen  such  that  the  possible 
error  bound  is  minimized.  This  research  is  the  first  to  discuss  this  issue. 


1 


2 


4.  Based  on  the  algorithm  for  the  single-item  problem,  this  research  also  provides 
a heuristic  algorithm  for  the  multi-item  dynamic  lot  size  model.  Preliminary 
computational  experience  shows  that  this  heuristic  algorithm  is  very  effective. 

The  outline  of  this  research  is  as  follows.  A literature  review  on  topics  of  dynamic 
lot  size  models  is  placed  in  Chapter  2.  In  this  chapter,  formulations  of  different  models 
are  stated,  and  algorithms  in  the  literature  for  them  are  discussed.  The  results  of 
this  research  are  in  Chapters  3 and  4.  Chapter  3 presents  the  results  related  to 
the  single-item  dynamic  lot  size  model,  including  an  algorithm  for  capacitated  lot 
size  models,  an  algorithm  for  minimizing  the  error  bound,  and  an  algorithm  for  the 
models  with  two  control  variables.  The  heuristic  algorithm  for  the  multi-item  models 
are  the  content  of  Chapter  4.  Conclusions  and  recommendations  for  future  research 
are  placed  in  the  last  chapter. 


CHAPTER  2 
LITERATURE  REVIEW 


2.1  Terminology  and  Notation  for  the  Lot  Size  Models 

Given  the  demands  for  a certain  number  of  periods,  cost  functions  of  produc- 
tion and  inventory,  and  technical  constraints  on  production  and  inventory  levels,  the 
dynamic  lot  size  model  seeks  to  determine  a feasible  lot  size  (or  production  level)  for 
each  period  with  the  total  cost  as  small  as  possible.  Here  “dynamic”  implies  that  the 
given  data  and  desired  decision  variables  vary  by  period  in  contrast  to  those  in  the 
traditional  EOQ  (economic  order  quantity)  model.  The  EOQ  model  assumes  that 
demands  and  costs  are  constant  over  time  implying  that  the  optimal  production  is 
also  at  a constant  level.  Thus  the  dynamic  lot  size  model,  being  more  realistic,  plays 
a very  important  role  in  production  planning  and  scheduling. 

2.1.1  Production  Planning  and  Scheduling 

“Production  planning  is  primarily  concerned  with  the  adaptation,  or  more  ex- 
actly the  tuning  of  the  firm’s  industrial  resources,  in  order  to  meet  demand  for  its 
final  products.”  (page  1,  Bensourssan,  Crouhy  and  Proth,  1983) 

“Production  schedules  are  simply  a means  of  translating  orders  into  a timetable 
for  manufacturing  activities.”  (page  17,  Buxey,  1989) 

Two  factors,  namely,  variability  of  demands  and  resources  and  capacity  of  re- 
sources, make  production  planning  and  scheduling  difficult.  In  general,  resources  of 
a company  include  capital  equipment,  raw  material,  labor,  and  time  available  to  fill 
demands,  and  all  of  these  resources  have  the  properties  of  limitation  and  variability. 
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Production  planning  can  be  divided  into  three  levels:  long-term,  medium-term 
and  short-term.  Long-term  production  planning  is  called  capacity  planning , which  is 
to  decide  the  type  and  the  number  of  machines  or  production  lines  based  on  long- 
term and  approximate  forecasting  of  demands.  In  the  medium-term,  it  is  expected 
that  the  forecast  demands  of  each  production  period  are  relatively  precise.  To  decide 
the  lot  size  of  each  period  is  the  main  topic  at  this  level.  Thus,  this  is  lot  size  plan- 
ning. The  third  level  is  machine  scheduling , which  concerns  how  to  assign  production 
jobs  to  proper  machines.  Often,  long-term  and  medium-term  planning  is  referred  as 
“production  planning”  and  short-term  planning  as  “production  scheduling.” 

The  concept  of  a planning  hierarchy  (Holstein,  1966)  can  also  be  used  to  illustrate 
production  planning  and  scheduling.  Typically,  production  planning  can  be  divided 
from  top  to  bottom  into  five  stages.  Each  stage,  which  is  limited  in  scope  by  its 
parent  stage,  considers  a plan  with  a shorter  time  span  and  more  detail  as  follows 
(Buxey,  1989). 

(1)  Capital  equipment  requirements  planning, 

(2)  Resources  (labor  and  material)  planning, 

(3)  Master  Production  Schedule  (MPS)  for  the  output  of  end  products, 

(4)  Materials  Requirement  Planning  (MRP)  for  the  extrapolation  of  MPS  to 
lower  levels  of  production, 

(5)  Machine  scheduling  for  actually  production. 

The  dynamic  lot  size  model  covers  stages  (2)  to  (4). 

2.1.2  Dynamic  Lot  Size  Models 

In  this  section,  we  will  discuss  the  terminology  and  some  properties  of  the  dy- 
namic lot  size  model.  We  will  consider  the  deterministic  version  of  this  model. 

Usually,  two  kinds  of  costs  are  considered.  The  production  cost  corresponds  to 
the  production  level,  and  the  inventory  cost  occurs  because  of  carrying  inventory.  The 
production  level  is  always  required  to  be  nonnegative.  The  inventory  level  may  be 
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nonnegative  or  negative.  The  latter  case  is  called  backlogging.  The  cost  functions  are 
always  assumed  to  be  nondecreasing.  Sometimes,  the  relationship  of  costs  between 
periods  is  considered.  When  a cost  is  said  to  be  nondecreasing  or  nonincreasing, 
it  means  that  the  unit  cost  is  constant  within  a period  but  is  nondecreasing  (or 
nonincreasing)  over  time  periods.  On  the  other  hand,  a concave  or  convex  cost  means 
it  is  concave  or  convex  over  production  levels  in  one  period. 

Capacitated  models.  One  reason  for  capacities  in  a problem  is  the  scarcity  of 
resources.  Another  reason  is  technical  limitations.  For  example,  every  machine  has 
its  design  speed.  It  cannot  be  operated  exceeding  this  speed  to  produce  more  product. 
It  is  important  to  consider  capacity  constraints  if  the  limitation  of  a certain  resource 
or  machine  may  be  reached.  Otherwise,  the  model  may  yield  an  infeasible  solution. 

Multi-item  models.  In  modern  industry,  almost  no  manufacturing  company 
produces  a single  product.  Besides,  a product  usually  is  assembled  from  many  indi- 
vidual parts.  A bill  of  materials  indicates  all  of  the  information  about  what  and  how 
many  parts  are  needed  to  assemble  a final  product.  Producing  these  parts  is  similar 
to  producing  multiple  items  of  products.  When  more  than  one  item  competes  for  a 
limited  resource,  it  becomes  a multi-item  problem. 

Multi-level  models.  The  name  “multi-level”  comes  from  the  bill  of  materials 
when  there  is  more  than  one  step  to  assemble  a final  product.  A multi-level  dynamic 
lot  size  problem  occurs  when  an  assembly  part  (partial  product)  is  also  a final  product 
or  can  be  carried  over  for  future  use.  For  example,  auto  parts  are  both  final  products 
and  assembly  parts  that  can  be  carried  over. 

Variations.  In  addition  to  the  above  models,  there  are  many  variations  of  cost 
functions  and  production  constraints.  When  a machine  or  production  line  is  set  up,  in 
addition  to  setup  cost,  there  may  exist  some  setup  time.  A learning  setup  means  that 
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setup  cost  or  time  is  a (decreasing)  function  of  the  number  of  setup.  A startup  cost 
means  that  a machine  will  cost  more  to  set  up  if  it  is  shut  down  for  several  periods.  In 
some  cases,  a reservation  cost  must  be  paid  to  keep  a machine  on  to  prevent  startup 
costs,  even  if  there  is  no  production.  In  the  multi-item  problems,  a change-over  cost 
may  occur  when  a machine  is  changed  from  one  product  to  produce  another  product 
in  the  same  period. 

Three  production  patterns  other  than  the  normal  fractional  production  have 
been  used  in  the  literature.  Full  batch  production  means  that  once  an  item  is  produced 
it  is  produced  in  an  amount  to  satisfy  the  exact  demands  of  certain  periods.  All- 
or-nothing  production  states  that  at  most  one  item  can  be  produced  each  period. 
Sometimes,  the  production  level  of  an  item  is  the  total  capacity  of  machines  assigned 
to  this  item.  Hence,  this  is  a problem  of  machine  assignment. 

2.1.3  Notation 

The  following  notation  is  used  in  the  remainder  of  this  chapter. 

T = time  horizon; 

/ = number  of  items; 
dt  = demand  in  period  t; 

Ct  = capacity  level  in  period  f; 
xt  = production  level  in  period  t; 

Dt  = cumulative  demand  up  to  period  f,  that  is,  Dt  = £*=1  dp, 

Ct  — cumulative  capacity  level  up  to  period  f,  that  is,  Ct  — SLi  cb 
Xt  = cumulative  production  level  up  to  period  t,  that  is,  Xt  = £j=1  xp 
yt  = setup  index  in  period  f; 

It  = inventory  level  in  period  t\ 

Pt(xt)  = production  cost  function  in  period  t; 

Ht{It)  = holding  (including  backlogging)  cost  function  in  period  t; 

Kt  = fixed  setup  cost  in  period  t; 
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Pt  = unit  production  cost  in  period  t ; 
ht  = unit  holding  cost  in  period  t; 


bt  = unit  backlogging  cost  in  period  t; 
dn  = demand  of  item  i in  period  f; 
xn  = production  level  of  item  i in  period  t\ 
yn  = setup  index  of  item  i in  period  t; 
lit  — inventory  level  of  item  i in  period  t; 

Pit  (xu)  = production  cost  function  of  item  i in  period  t ; 

Hit(Itt)  = holding  (including  backlogging)  cost  function  of  item  i in  period  f; 
Kn  — fixed  setup  cost  of  item  i in  period  f; 

Pit  = unit  production  cost  of  item  i in  period  t; 
ha  = unit  holding  cost  of  item  i in  period  t; 

= fixed  (setup)  resource  consumption  of  item  i in  period  t; 
fat  = unit  resource  consumption  rate  of  item  i in  period  t; 
rij  = number  of  item  i required  to  produce  one  unit  of  item  j. 


2.2  Single-Item  Dynamic  Lot  Size  Models 


In  the  following  models,  there  is  no  initial  and  ending  inventory: 


Io  = 0,  It  — 0. 


(2.1) 


lliere  is  no  loss  of  generality  in  this  assumption.  If  there  is  initial  inventory,  reduce 
the  demand  of  the  first  period  required  (or  those  of  the  first  several  periods)  by  that 
amount.  If  there  is  ending  inventory,  increase  the  demand  of  the  last  period(s)  by 
that  amount. 


Define  Pt(xt)  and  Ht(It)  to  be  the  production  cost  and  holding  cost  functions  in 
period  t.  Usually  Pt(xt ) is  defined  as: 


0 if  xt  = 0 

I<t  + ptxt  if  xt  > 0 


(2.2) 
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where  K%  is  the  setup  cost  and  pt  is  the  unit  production  cost  for  any  1 < i < T. 
Holding  cost  is  defined  as 


-btIt  if  7*  < 0 
htIt  if  It  >0 


where  bt  is  the  unit  backlogging  cost  and  ht  is  the  unit  normal  holding  cost  for  any 
1 < t < T.  Note  that  if  backlogging  is  not  allowed,  Ht(It)  is  only  defined  over  > 0 
for  1 < i < T. 


2.2.1  Uncapacitated  Models 


Wagner  and  Whitin  (1958)  first  studied  the  dynamic  lot  size  model  (DLSM). 
The  DLSM  can  be  formulated  as 


T 

min  Y^KtUt  + Ptxt  + htIt 

t= 1 


such  that 


It- 1 + Xf  — df  + V t 

(2.4) 

o 

II 

-£r 

ii 

(2.5) 

T 

xt<(52dk)yt,  Vi, 

(2.6) 

k=t 

Xt,  It  > 0,  Vi 

(2.7) 

yt  € {0, 1},  Vi 

(2.8) 

When  pt  = 0 for  all  t,  it  is  called  the  Wagner- Whitin  model 

Zangwill  (1966)  extended  the  Wagner- Whitin  model  to  include  backlogging. 
Notice  that  the  only  difference  between  problems  with  and  without  backlogging  is 
whether  the  model  includes  the  nonnegative  constraints  of  inventory  levels. 

In  the  literature,  there  are  two  basic  types  of  solution  approaches  to  the  DLSM. 
They  are  dynamic  programming  and  linear  programming  with  facet  cuts.  Below  these 
will  be  briefly  summarized. 
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Dynamic  programming.  Wagner  and  Whitin  (1958)  gave  the  following  property 
of  optimality  in  their  pioneering  paper. 

Lemma  1 (property  of  optimality)  There  is  an  optimal  policy  for  DLSM  such 
that  xtIt_ i = 0 for  any  t.  That  is,  in  any  period  where  there  is  production  the 
incoming  inventory  must  be  zero. 

The  property  of  optimality  implies  the  lemma  below: 

Lemma  2 If  {xt)  is  an  optimal  policy  for  DLSM,  then  every  xt  = 0 or  xt  = di 
for  some  k between  t and  T. 

Based  on  the  above  lemma,  Wagner  and  Whitin  provided  a dynamic  programming 
solution  with  complexity  0(T2).  Let  F(t)  denote  the  minimal  cost  to  schedule  pro- 
duction for  periods  1 through  t with  It  = 0.  With  F( 0)  = 0, 

F(t)  — min  {F(s  — 1)  + Ks  -f  H(s,t)} for  t>  1, 

where 

t t t 

H(s,  t ) = hs  dk  + hs. |-i  dk  + . . . + ht+i  dk,  for  s < t. 

k=s+l  k=s-\-  2 k=t 

Veinott  (1968)  extended  the  argument  of  Wagner  and  Whitin  by  showing  that 
the  lemmas  hold  for  any  concave  cost  function.  Since  then,  this  problem  has  been 
extensively  studied  by  several  other  authors  including  Love  (1973),  Swoveland  (1975), 
Johnson  and  Montgomery  (1974),  Florian  et  al.  (1980)  and  Silver  and  Peterson 
(1985). 

Recently,  Wagelmans  et  al.  (1992)  explored  a geometric  property  and  provided  a 
backward  dynamic  programming  solution  with  computational  complexity  O(TlogT). 
This  was  the  first  paper  to  improve  the  complexity  bound  since  Wagner  and  Whitin. 
Independently,  Fedurgruen  and  Tzur  (1991)  and  Aggarwal  and  Park  (1993)  provided 
algorithms  that  also  work  in  0(T  log  T).  In  particular,  Fedurgruen  and  Tzur  used 
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a forward  dynamic  programming  approach,  and  the  Aggarwal  and  Park  algorithm 
took  advantage  of  matrix  computation  instead  of  a dynamic  program.  Motivated  by 
Wagelmans  et  al.  (1992),  Chen  and  Lee  (1991)  gave  a forward  dynamic  program  with 
the  same  computational  complexity.  Fedurgruen  and  Tzur  (1991)  also  showed  that 
the  DLSM  with  backlogging  can  be  solved  in  0(T log  T). 

Linear  programming  with  facet  cuts.  Attempting  to  use  linear  programming 
to  solve  the  DLSM  is  natural,  since  it  has  already  been  formulated  as  a mathematical 
program. 

A facet  cut  is  an  inequality  constraint  that  forms  a facet  of  the  convex  hull  of 
all  feasible  integer  solutions.  If  all  the  facet  cuts  are  found,  then  the  solution  of  the 
LP  relaxation  will  be  exactly  the  solution  of  the  mixed  integer  program. 

Let  dij  — J2t=i  dt-  Let  L = {1, . . . , /}  for  any  1 < / < T,  and  let  S C L.  Barany 
et  al.  (1984a)  showed  that  the  following  are  valid  inequalities  for  the  previous  model. 

^2  xi  + H duyi  > d,u  (2.9) 

ies  i£L\S 

Moreover,  they  also  showed  that  constraints  (2.4)  - (2.7)  plus  (2.9)  and  0 < yt  < 1 
for  all  t form  the  integer  convex  hull  of  constraints  (2.4)  - (2.8). 

Note  that  this  requires  adding  an  exponential  number  of  constraints  to  the  previ- 
ous model.  Therefore,  the  facet  cut  approach  does  not  provide  an  efficient  algorithm 
for  the  Wagner- Whitin  model.  However,  research  on  the  multi-item  problem  has 
shown  that  inequalities  similar  to  (2.9)  are  useful  on  such  problems.  Papers  applying 
this  idea  can  be  found  in  Barany  et  al.  (1984b),  Leung  et  al.  (1989)  and  Pochet  and 
Wolsey  (1991). 

2.2.2  Planning  Horizon  and  Error  Bound 

The  number  of  periods  in  the  dynamic  lot  size  model  is  called  the  time  horizon. 
There  are  two  kinds  of  time  horizons.  The  number  of  periods  for  which  planning 
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must  occur  is  called  the  natural  horizon ; by  contrast,  the  number  of  periods  for  which 
the  decision  maker  has  reliable  data  about  demands  and  costs  is  call  the  data  horizon 
(Denardo  and  Lee  (1991)).  Obviously,  the  lot  size  decision  is  based  on  the  data 
horizon.  One  question  arises:  “How  long  should  the  data  horizon  be  to  make  the 
right  decision  for  lot  sizes  of  the  first  several  periods?”  The  time  horizon  that  is  long 
enough  to  make  such  a decision  is  called  a forecast  horizon  or  a planning  horizon.  A 
perfect  horizon  is  the  smallest  forecast  horizon  needed. 

In  the  dynamic  lot  size  model,  a best  solution  for  the  natural  horizon  is  sought. 
Since  the  DLSM  is  a deterministic  model,  the  optimal  solution  can  be  found  if  all 
the  information  (perfect  information)  is  available  for  the  natural  horizon.  In  the  real 
world,  however,  data  are  available  only  for  the  data  horizon,  which  is  usually  far 
less  than  the  natural  horizon.  Unfortunately,  the  first  decision  must  be  made  under 
this  imperfect  information.  Here  the  first  decision  is  the  production  plan  for  the  first 
several  periods.  Once  the  first  decision  has  been  executed,  it  cannot  be  changed  even 
if  it  is  later  learned  this  not  a best  decision  because  of  new  information.  Therefore, 
the  concern  is  whether  the  information  at  hand  is  enough  to  make  the  first  decision 
and  how  much  error  may  occur  because  the  first  decision  is  made  under  this  imperfect 
information.  The  maximal  error  that  occurs  because  the  first  decision  is  made  under 
imperfect  information  is  called  an  error  bound. 

Literature  in  this  field  addresses  four  things:  the  planning  horizon,  the  per- 
fect horizon,  and  existence  conditions  for  the  planning  horizon  and  the  error  bound. 
Wagner  and  Whitin  (1958)  discussed  the  planning  horizon  in  their  paper.  Later  it 
was  extended  by  Zabel  (1964),  Kunreuther  and  Morton  (1973,  1974),  Blackburn  and 
Kunreuther  (1974),  Lundin  and  Morton  (1975)  and  Bean  et  al.  (1987).  Chand  and 
Morton  (1979,  1986)  provided  algorithms  for  the  perfect  horizon.  Existence  condi- 
tions for  the  planning  horizon  are  discussed  in  the  papers  by  Eppen  et  al.  (1969), 
Bean  and  Smith  (1984,  1990)  and  Chand  et  al.  (1990).  Lee  and  Denardo  (1986)  and 
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Denardo  and  Lee  (1991)  provided  an  0{T 2)  method  to  determine  the  error  bound 
for  the  problem  with  nonincreasing  production  costs.  Moreover,  this  algorithm  needs 
data  beyond  period  T.  Federgruen  and  Tzur  (1991)  gave  an  O(riogT)  algorithm  to 
find  an  upper  bound  for  the  error  bound. 


2.2.3  Capacitated  Models 

The  single-item  capacitated  dynamic  lot  size  model  (CLSM)  is  similar  to  the 
DLSM,  but  there  is  some  limitation  of  a scare  resource,  such  as  production  capacity 
or  warehouse  space.  Let  ct  be  the  limitation  of  production  level  in  period  t.  The 
CLSM  formulation  is  just: 


T 

minJ^  KtUt  + Pt%t  + hflt 
t= i 


such  that 


It- 1 + %t  — dt 

+ /t,  Vt 

(2.10) 

ii 

ii 

o 

(2.11) 

Xt  < CtVt, 

Vt 

(2.12) 

Xti  It  S 0, 

vt 

(2.13) 

Vt  s {o?i}, 

Vt 

(2.14) 

Similarly,  let  ct  be  the  inventory  capacity  in  period  t.  Then  the  CLSM  with  inventory 
capacity  is  as  follows. 

T 

min  YjKtVt  + Ptxt  + htIt 

t= i 

such  that 


It- i + Xt  — dt  T 7t,  V t 

(2.15) 

Iq  = It  = 0, 

(2.16) 

T 

xt  < (£dk)yt,  Vt 

(2.17) 

k=t 
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h <cu  V t 

(2.18) 

Xt  < Ctyt,  V t 

(2.19) 

xt,  It>  0,  V t 

(2.20) 

yt  E {0, 1},  V t 

(2.21) 

The  computational  complexity  of  the  CLSM  has  been  investigated  by  Florian  et 
al.  (1980)  and  later  by  Bitran  and  Yanasse  (1982).  They  showed  that  the  CLSM  is 
NP-hard  even  in  many  special  cases.  Therefore,  researchers  have  tried  to  define  special 
cases  for  which  polynomial  algorithms  can  be  developed.  The  notation  a/f3/^/8, 
which  was  introduced  by  Bitran  and  Yanasse  (1982)  to  describe  families  of  the  CLSM, 
where  a , /?,  7 and  8 represent  setup  cost,  holding  cost,  production  cost  and  capacity 
type,  respectively.  The  value  of  these  parameters  may  be  G , C,  ND , NI  and  Z for 
general  pattern,  constant,  nondecreasing,  nonincreasing  and  zero.  In  this  notation, 
the  production  costs  are  linear.  Florian  and  Klein  (1971)  gave  an  0(T4)  algorithm 
for  G/G/G/C . Jagannathan  and  Rao  (1973)  extended  this  result  to  general  period 
cost  functions.  Swoveland  (1975)  adapted  Florian  and  Klein’s  algorithm  for  piecewise 
concave  period  cost  functions.  Love  (1973)  showed  that  the  same  model  with  constant 
inventory  capacity  instead  of  production  capacity  is  0(T3).  Bitran  and  Yanasse 
showed  that  the  time  complexity  of  NI/G/NI/ND , NI/G/NI/C , C/Z/C/G  and 
ND/Z/N D /NI  are  0(T4),  0(T3),  0(T  logT)  and  0(T),  respectively.  Chung  and 
Lin  (1988)  improved  the  time  complexity  of  NI/G/NI/ND  to  0(T2). 

For  the  general  CLSM,  Florian  et  al.  (1980)  suggested  a pseudo-polynomial 
dynamic  programming  formulation  with  complexity  of  0(DTCT),  where  DT  is  the 
total  demand  and  Cy  is  the  total  capacity.  Baker  et  al.  (1978)  provided  a tree- 
search  solution  algorithm  for  the  G/G/C/G  cases.  For  problems  with  G/G/G/G , 
there  are  two  papers,  by  Lambrecht  and  Van  der  Eecken  (1978)  and  by  Kirca  (1990). 
The  former  implemented  the  optimal  property  from  Florian  and  Klein  (1971)  and 
designed  a dynamic  programming  solution  procedure.  The  latter  refined  the  idea  of 
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pseudo-polynomial  dynamic  programming.  Kirca’s  algorithm  seems  to  be  the  better 
of  the  two,  however,  it  needs  328  seconds  (on  CORDATA-ATP-8-Q20)  to  solve  an 
18-period  problem  for  the  worst  case  out  of  240  tests.  The  computing  time  of  this 
algorithm  increases  at  a rate  of  0(T4).  Chung  et  al.  (1990)  provided  an  algorithm 
that  combines  dynamic  programming  and  branch-and-bound.  This  algorithm  only 
works  on  the  problem  with  nondecreasing  production  costs. 

2.2.4  Other  Models 

Startup  and  reservation  costs.  Karmarkar  et  al.  (1987)  considered  the  CLSM 
with  startup  cost  (CLSM/S).  They  split  conventional  setup  cost  in  the  CLSM  into 
startup  cost  and  reservation  cost.  A startup  cost  occurs  when  a machine  is  on  in  one 
period  and  off  in  the  previous  period.  A reservation  cost  occurs  once  a machine  is  on. 
Let  St  be  the  startup  cost  and  Kt  the  reservation  cost.  This  model  can  be  formulated 
as  follows. 

T 

min  StZt  + Ktyt  + ptxt  + htIt 


t= 1 


such  that 


It- 1 + V t 

(2.22) 

II 

ii 

o 

(2.23) 

Xt  < CtVu  v t 

(2.24) 

zt>yt-yt- 1,  V t 

(2.25) 

xt,It  >0,  V t 

(2.26) 

yt,zt  e {0,1},  Vf 

(2.27) 

This  model  is  similar  to  the  CLSM  if  the  reservation  cost  is  considered  as  the  setup 
cost  in  the  latter  model.  Therefore,  the  only  difference  is  that  one  more  cost  (startup 


cost)  is  considered  in  the  CLSM/S. 
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Karmarkar  et  al.  used  dynamic  programming  for  the  uncapacitated  case  and  pro- 
vided a branch-and-bound  algorithm  with  Lagrangean  relaxation  for  the  capacitated 
case.  The  empirical  performance  of  this  algorithm  is  0(T4).  A dynamic  program, 
which  is  in  the  complexity  of  0(T2),  was  also  provided  in  the  same  paper  for  the 
uncapacitated  problem.  Wolsey  (1989)  applied  a facet  cut  approach  to  this  problem. 


Learning  in  setups.  Chand  and  Sethi  (1990)  assumed  that  setup  costs  in  the 
DLSM  are  not  independent  of  each  other,  but  a function  of  the  number  of  setups. 

That  is,  Kt  = K (St)  where  St  = ^Li  <HX«)-  The  dynamic  lot  size  model  with  learning 
in  setups  (DLSM/L)  is  as  follows. 

T 

£*(«.)  »<  + PtXt  + htIt 
t= i 


such  that 


It- 1 + xt  — dt 

+ h,  yt 

(2.28) 

S'* 

II 

•s" 

II 

o 

VJ 

(2.29) 

t,  i L ^ 0, 

yt 

(2.30) 

yt  € {o,i}, 

yt 

(2.31) 

where  6t  = Xw=i  Vi  and  K(S)  is  a nonincreasing  function. 

Karwan  et  al.  (1988)  and  Chand  (1989)  discussed  the  learning  effect  for  the 
cases  of  infinite  horizon  and  constant  demand  rate.  Chand  and  Sethi  (1990)  provided 
a dynamic  programming  algorithm  for  the  DLSM/L  with  complexity  of  0(T3). 


Multiple  sources.  In  the  previous  discussion,  it  was  assumed  that  demand  can 
only  be  satisfied  by  a single  source.  This  is  not  always  the  case.  For  example,  rather 
than  producing  it  may  be  possible  to  buy  from  an  outside  vendor.  In  another  case,  a 
big  company  may  have  many  factories  in  different  locations  with  different  production 
costs.  Lee  and  Zipkin  (1989)  considered  a two-source  problem  with  the  decision  of 
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producing  or  buying.  Let  zt  be  the  amount  to  buy  with  the  unit  cost  qt  in  period  t. 
Assume  that  there  is  no  capacity  limitation  when  buying.  The  capacitated  dynamic 
lot  size  model  with  multi-sourcing  (CLSM/S)  is  formulated  as  follows. 

T 

min  Ktyt  + pt%t  + Qt^t  + hflt 

(=i 


such  that 


It- 1 + %t  + Zt  — dt  + /*,  Vi 

(2.32) 

II 

II 

O 

(2.33) 

Xt  < Ctyt,  V t 

(2.34) 

Xti  It  ^ o,  V t 

(2.35) 

yt  € {0, 1},  V t 

(2.36) 

Erenguc  and  Tufekci  (1987)  used  the  Lagrangean  relaxation  technique  with 
branch  and  bound.  Problems  with  10  periods  were  solved  in  about  1 second  on 
the  IBM3081/D  (1.7  Mflops  per  second).  Lee  and  Zipkin  (1989)  provided  an  0(T 5) 
algorithm  for  the  case  of  constant  production  capacity. 

2.3  Multi-item  Dynamic  Lot  Size  Models 

The  multi-item  lot  size  problem  occurs  when  there  is  more  than  one  production 
item  competing  for  a scarce  resource.  For  example,  several  products  are  produced 
from  the  same  machine  or  use  the  same  warehouse,  where  the  machine  has  its  time 
capacity  and  the  warehouse  has  its  space  limitation. 

2.3.1  The  Model  with  a Common  Resource  Capacity 

In  this  subsection,  we  will  discuss  the  multi-item  dynamic  lot  size  model  with 
one  common  capacity  constraint  (MLSM).  The  notation  used  here  is  similar  to  that 
of  the  last  section  but  with  an  additional  subscript  i to  indicate  item  i.  For  example, 
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xu  means  the  production  level  of  item  i in  period  t.  Let  the  total  number  of  items  be 
I.  The  MLSM  is  formulated  below. 

I T 

Kityn  + Pit x a + hula 

2 — 1 t= 1 

such  that 


Ii,t— 1 "f"  x it  d{t  “f"  lit 5 V 

T 

(2.37) 

1 

fcXit  - v t 

i= l 

(2.38) 

lio  = IiT  — 0,  V i 

T 

(2.39) 

xit  ^ (^  ^ djj )y^7 
j—t 

(2.40) 

'll iti  ^ 0,  V If  t 

(2.41) 

Vit  € {0, 1),  V i,  t 

(2.42) 

where  fa  is  the  resource  consumption  rate  of  item  i. 

There  are  / x T zero-one  variables  in  the  above  model.  For  a half-year  plan 
with  weekly  lots  and  100  items,  the  total  number  of  integer  variables  will  be  about 
2500,  making  it  impractical  to  solve  this  model  directly  by  a mixed  integer  linear 
program.  In  the  rest  of  this  section,  we  will  discuss  common  techniques  proposed  in 
the  literature  for  this  model. 

Lagrangean  relaxation  and  branch-and-bound.  Usually,  constraint  (2.38)  is  re- 
laxed in  the  MLSM.  Let  7r(  be  the  corresponding  multipliers.  The  Lagrangean  dual 
is  as  follows. 

it  t 

L(tt)  = min^  £ (KitVn  + ( Pit  + fat)xit  + hitIit)  + ^ ctTTt 

i= 1 t= 1 t= 1 

such  that 

Ii,t-i  + xit  — dn  + In ->  V i,t  (2.43) 

lio  = IiT  = 0,  V i (2.44) 
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T 

XH  < E^,  V i,t  (2.45) 

j=t 

Xu,  lit  > o,  Vi,f  (2.46) 

3/»t  € {0,1},  V i,t  (2-47) 

Then  T(7r)  can  be  separated  into  / single-item  problems  with  a particular  7 r given. 
The  dual  problem  is  to  maximize  T(7r)  with  7r  nonnegative.  It  is  well  known  that 
L(tt)  is  a concave  function  (Bazaraa  and  Shetty,  1979). 

Normally,  three  issues  are  of  concern  when  applying  this  relaxation  scheme: 

(1)  How  to  create  the  dual  problem,  or  equivalently,  which  constraints  to  relax. 

(2)  How  to  solve  the  dual  problem.  This  issue  includes  the  algorithm  for  the 
relaxed  (single-item)  subproblems  and  the  method  to  update  the  multipliers. 

(3)  How  to  recover  a (primal)  feasible  solution  or  what  is  the  branching  rule  if 
a branch-and-bound  process  is  used  to  find  an  optimal  solution. 

Most  papers  in  the  literature  choose  some  version  of  the  subgradient  method 
to  update  the  multipliers.  Thizy  and  Van  Wassenhove  (1985)  and  Gelders  et  al. 
(1986)  used  a subgradient  method  to  update  multipliers  and  solved  a transportation 
problem  to  get  an  upper  bound.  The  convergence  rate  was  shown  to  be  very  slow 
near  an  optimal  solution.  For  example,  one  of  the  8 x 8 (I  x T ) test  problems  in  the 
second  paper  above  needed  69  seconds  on  the  IBM  4341  with  1.7  Mflops  per  second 
for  an  optimal  solution  but  only  0.66  seconds  for  an  approximate  solution  with  1% 
error.  Hence,  research  is  often  directed  at  finding  an  approximate  solution.  Four 
of  Thizy  and  Van  Wassenhove’s  8x8  problems,  which  need  CPU  time  for  optimal 
solutions  ranging  from  4290  seconds  to  69  seconds  (on  the  IBM  4341),  have  become 
benchmarks  for  many  proposed  heuristic  methods. 

LP  relaxation  and  facet  cuts.  Manne  (1958)  was  first  to  use  linear  program- 
ming on  dynamic  lot  size  models.  Suppose  that  5,-  is  the  set  of  possible  plans  for 
item  i where  Si  is  based  on  the  Wagner- Whitin  property  of  optimality.  Since  in  each 
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period  the  amount  produced  is  equal  to  the  demand  of  the  next  several  periods,  there 

are  a total  of  2T  1 of  these  plans  for  each  item.  Let  pij  be  the  corresponding  cost  of 

plan  j for  item  i and  similarly  C{jt  be  the  resource  needed  at  period  t.  Introduce  the 
variable 


q f 1 if  schedule  j is  used  for  item  z, 

u [ 0 otherwise. 

Then  the  MLSM  can  be  formulated  as 

/ 

minL  PiAi 

*=1  i€5i 

such  that 


(2.48) 


£%  = 1,  Vi 

(2.49) 

jeSi 

T 

1 

J2  CiP0ij  -c*i  v i 

(2.50) 

i= i jeSi 

%>  0,  V i,j 

(2.51) 

Later,  Dzielinski  and  Gomory  (1965)  implemented  Dantzig- Wolfe  decomposition  to 
avoid  enumerating  all  possible  plans.  Lasdon  and  Terjung  (1971)  used  a similar 
column  generation  technique  to  this  model.  Bahl  (1983)  gave  a column  generation 

based  heuristic  algorithm.  Aardal  and  Larsson  (1990)  applied  Benders  decomposition 
to  this  model. 

The  key  drawback  of  this  approach  is  that  only  the  solutions  that  satisfy  the 
optimality  property  for  the  single-item  model  are  obtained.  However,  this  property 
may  not  hold  in  the  solution  for  the  MLSM. 

Eppen  and  Martin  (1987)  applied  a variable  redefinition  concept  to  MLSM.  Vari- 
able redefinition,  in  short,  is  to  describe  the  equivalent  network  flow  problem,  usually 
a shortest  path  problem,  instead  of  the  original  DLSM.  The  gap  of  LP  relaxation 
is  reduced  because  of  a tighter  formulation.  The  disadvantage  is  that  it  needs  IT 2 
decision  variables.  This  algorithm  solved  problems  with  100  items  and  10  periods  in 
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about  30  minutes  on  a Vax  11/750  (0.12  Mflops  per  second)  with  0.3%  error  compared 
to  best  known  solutions. 

Unlike  the  above  approach,  which  required  many  variables,  Barany  et  al. 
(1984a, b),  Pochet  and  Wolsey  (1988,  1991)  and  Leung  et  al.  (1989)  attempted  to 
add  facet  cuts  to  reduce  the  duality  gap.  Pochet  and  Wolsey’s  algorithm  took  about 
250  seconds  on  a Data  General  MV8000  (0.078  Mflops  per  second)  to  solve  the  linear 
programs  with  cuts  for  problems  with  5 items  and  12  periods. 

Heuristic  approaches.  The  MLSM  is  known  to  be  NP-hard  since  the  single- 
item case  is  NP-hard  already  (see  Florian  et  al.  (1980)  and  Bitran  and  Yanasse 
(1982)).  Therefore,  many  heuristic  algorithms  were  designed  for  this  model.  We 
cite  Eisenhut  (1975),  Lambrecht  and  Vanderveken  (1979),  Van  Nunen  and  Wessels 
(1978),  Dixon  and  Silver  (1981),  Dogramaci  et  al.  (1981),  Kami  and  Roll  (1982),  Maes 
and  Van  Wassenhove  (1986a, b)  and  Thizy  and  Van  Wassenhove  (1985).  An  excellent 
review  was  done  by  Maes  and  Van  Wassenhove  (1988).  Most  of  the  heuristic  methods 
try  to  allocate  resource  profitably  period  by  period  and/or  to  improve  a feasible  plan 
at  hand  by  some  “common  sense.”  Thizy  and  Wassenhove  algorithm  is  the  only 
one  that  used  an  optimization  model.  They  solved  a Lagrangean  relaxation  problem 
for  the  lower  bound  and  solved  a transportation  problem  at  each  iteration  to  find  a 
feasible  solution  as  the  upper  bound.  This  algorithm  gave  the  best  solution  (least 
erroi  from  optimal  solution)  but  needed  more  computational  effort  than  other  ones. 
For  example,  this  algorithm  found  optimal  solutions  for  9 out  of  their  14  test  problems 

in  50  iterations,  and  its  average  CPU  time  is  200  times  more  than  that  of  Dixon  and 
Silver’s  algorithm. 

2.3.2  Models  with  Setup  Time 

In  the  real  world,  to  set  up  a machine  requires  both  a fixed  charge  cost  and  a 
fixed  amount  of  resource  consumption  (setup  time).  Let  ccj  be  the  setup  time  required 
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for  item  i.  The  multi-item  dynamic  lot  size  model  with  setup  time  (MLSM/T)  can 
be  formulated  as  follows. 

I T 

min  ^ ^ KitVit  + Pit x it  + hit  lit 

2—1  t= 1 

such  that 


I{,t- 1 + xit  — dit  + lit,  V i,t 

T 

(2.52) 

1 

X ( + fan)  <ct,  V t 

2 = 1 

(2.53) 

* 

> 

o 

II 

ii 

1 o 
. ^ 

(2.54) 

xit  — ( ^ y djj  ^)llit , v l^t 

j—t 

(2.55) 

J 

x it')  lit  ^9?  V 

(2.56) 

Hit  G {0, 1},  V id 

(2.57) 

Newson  (1975a, b),  Billington  et  al.  (1986)  and  Trigeiro  (1989)  have  heuristic 
algorithms  for  this  model.  Moreover,  Maes  et  al.  (1991)  showed  that  finding  a feasible 
solution  is  NP-complete  for  problems  including  setup  times.  So,  no  feasible  solution 
is  guaranteed  in  the  existing  methods. 

2.3.3  Multi-level  Models 

A multi-level  lot  size  problem  occurs  when  an  item  serves  as  both  a final  product 
and  an  assembly  part.  Let  Si  be  the  set  of  successors  of  item  i and  r2-j  be  the  number 
of  units  of  item  i required  to  produce  one  unit  of  j . The  multi-level  dynamic  lot  size 
model  (MLSM/L)  is  as  follows. 

I T 

mm  ^ ^ I^itPit  + Pit'Eit  "T  hit  I it 

i=i  t= i 

such  that 


Ii,t-1  + Xit  - lit  = dit  + X 

jes, 

T 

V i,t 

(2.58) 

1 

X]  {otiVit  + fait)  < C<,  V t 

i= 1 

(2.59) 
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Iio  = IiT  = 0,  Vi  (2.60) 

T 

it  — y 1 ) !Ji t ' (2.61) 

j=t 

xiti  lit  >0,  V t,  t (2.62) 

Vit  e {0,1},  V*,<  (2.63) 


The  MLSM/L  without  constraint  (2.59)  is  called  the  uncapacitated  multi-level 
lot  size  problem.  Afentakis  et  al.  (1984)  and  Afentakis  and  Gavish  (1986)  refor- 
mulated the  constraint  matrix  of  this  problem  to  establish  a problem  with  a block- 
diagonal  structure  linked  by  a set  of  coupling  constraints.  The  tested  problems  had 
one  end  item,  12  periods  and  40  stages.  Kuik  and  Salomon  (1990)  used  a simulated 
annealing  heuristic  for  this  model. 

Maes  et  al.  (1991)  provided  an  algorithm  based  on  rounding  an  LP  solution 
for  the  problem  without  setup  times.  Problems  up  to  3 end  items,  10  periods  and 
3 levels  were  tested.  A heuristic  algorithm  solved  these  problems  in  an  average  55.9 
seconds  on  a Vax  11/785  (0.18  Mflops  per  second)  with  only  1%  error  compared  to 
the  optimal  solutions. 

2.3.4  Other  Models 

Since  the  MLSM  is  difficult  to  solve,  some  authors  have  attempted  to  specify 
the  pattern  of  the  solution  in  order  to  make  the  problem  easy. 

The  first  approach  is  called  full  batch  production.  Full  batch  production  assumes 
that  once  an  item  is  set  up,  it  produces  the  amount  satisfying  the  exact  demands  up 
to  a certain  period.  Manne’s  (1958)  model  is  in  this  class. 

Fleischmann  (1990)  made  an  assumption  that  at  most  one  item  can  be  produced 

at  each  period  and  the  lot  size  must  be  the  full  capacity  of  that  period.  This  is  called 
an  all-or-nothing  production. 

The  third  approach  is  called  machine  assignment.  The  decision  variables  are 
defined  to  be  the  number  of  machines  assigned  to  each  item  at  each  period.  The  lot 
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size  of  an  item  is  the  summation  of  full  capacity  of  all  the  machines  assigned  to  that 
item.  Bitran  and  Matsuo  (1986b)  and  Eppen  and  Martin  (1987)  used  this  assumption 
in  their  research. 


CHAPTER  3 

SINGLE-ITEM  DYNAMIC  LOT  SIZE  MODELS 

3.1  New  Dynamic  Programming  Algorithm 

The  single-item  dynamic  lot  size  model  (DLSM)  can  be  described  as  follows.  For 
a finite  time  horizon  T , demand  for  a single  commodity  of  each  production  period 
is  known.  These  demands  may  be  satisfied  by  production  in  the  current  or  future 
periods,  by  inventories  from  previous  periods,  or  by  purchases  from  outside  suppli- 
ers. There  may  be  certain  kinds  of  capacity  restraints  on  production,  inventories  or 
purchases.  We  consider  costs  that  are  incurred  by  producing  products,  purchasing 
products,  or  holding  inventories.  Here  the  holding  costs  may  occur  because  of  either 
positive  or  negative  inventory,  with  the  latter  representing  backlogging.  The  DLSM 
seeks  a feasible  production  plan  with  the  total  cost  as  small  as  possible. 

In  this  study,  we  propose  a dynamic  programming  algorithm  for  dynamic  lot 
size  models  where  the  production  and  inventory  cost  functions  are  only  assumed  to 
be  piecewise  linear.  In  particular,  there  are  no  assumptions  of  convexity,  concavity 
or  monotonicity.  Thus  the  algorithm  addresses  most  variants  of  the  DLSM  appearing 
in  the  literature. 

In  addition  to  the  literature  mentioned  in  the  previous  chapter,  certain  variations 
of  the  DLSM  have  been  treated  by  other  authors.  For  example,  Lippman  (1969) 
considered  the  DLSM  with  the  decision  variables  being  the  number  of  trucks  to  deliver 
orders.  This  idea  was  extended  in  the  works  of  Sethi  and  Chand  (1981)  and  Chand 
and  Sethi  (1983)  for  the  capacitated  DLSM.  Multi-source  versions  of  the  capacitated 
DLSM  are  treated  by  Lee  and  Zipkin  (1989)  who  allowed  make-or-buy  decisions  and 
gave  an  0(T5)  algorithm  for  the  case  of  equal  production  capacities.  And  Erenguc 
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and  Tufekci  (1987)  and  Erenguc  and  Aksoy  (1990)  have  implemented  branch-and- 
bound  algorithms  for  the  capacitated  DLSM  with  piecewise  linear  production  cost 
functions. 

Of  the  above  papers,  only  a few  report  computational  results  with  cpu  times. 
Po  indicate  the  state-of-the-art,  we  cite  those  and  quote  the  cpu  times  of  the  most 
difficult  problem  in  each:  Kirca  (1990)  solved  a G/G/G/G  problem  with  18  periods 
in  328  seconds  on  a CORDATA-ATP-8-Q20;  Chung  et  al.  (1990)  solved  one  of  their 
G/C/C/G  problems  with  96  periods  in  103.47  seconds  on  a VAX  11/750;  Erenguc 
and  Tufekci  (1987)  solved  a 10-period  capacitated  DLSM  with  4 segments  in  the 
production  cost  functions  in  4.65  seconds  on  an  IBM3081/D32;  and  Erenguc  and 
Aksoy  (1990)  solved  a similar  problem  with  12  periods  and  2 segments  in  0.284  seconds 
on  an  IBM3090/400. 

In  section  4 we  give  extensive  computational  results,  including  Vax  8600  cpu 
times,  for  the  method  of  this  paper.  From  these,  it  can  be  estimated  that  our  cpu 
times  for  the  above  four  cases  would  be  less  than  0.08,  2.91,  0.34  and  0.09  seconds, 
respectively.  The  Vax  8600  is  rated  at  0.48  megaflops  by  Dongarra(1989)  while  the 
Vax  11/750,  IBM3081  and  IBM3090  have  ratings  of  0.12,  2.1  and  16,  respectively. 
(The  CORDATA-ATP-8-Q20  is  not  rated.)  Thus,  even  allowing  for  the  different 
computers,  we  feel  this  study  confirms  the  effectiveness  of  the  method  relative  to 
those  proposed  thus  far. 

Our  algorithm  is  motivated  by  the  approach  of  Wagelmans  et  al.  (1992)  who 
have  recently  improved  the  complexity  of  the  original  Wagner- Whitin  {G / G / G / oo) 
model  to  O(TlogT').  (Independently,  Federgruen  and  Tzur  [1991]  and  Aggarwal  and 
Park  [1990]  have  obtained  the  same  complexity  with  different  approaches.)  It  also 
represents  a significant  generalization  and  extension  of  our  work  on  the  capacitated 
DLSM  in  Chen  et  al.  (1992). 
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The  algorithm  results  from  a continuous  version  of  the  standard  dynamic  pro- 
gramming model.  The  state  variable,  that  is,  the  cumulative  production  level,  is 
treated  as  continuous.  Under  the  assumption  that  the  cost  functions  (production  and 
inventory)  are  piecewise  linear,  we  develop  an  efficient  method  to  update  the  optimal 
value  function.  While  the  computational  effort  depends  on  the  number  of  segments  in 
the  optimal  value  function,  which  can  increase  exponentially  with  T,  extensive  testing 
has  shown  the  method  to  be  very  effective  over  a wide  range  of  problem  parameters. 

3.1.1  Notation.  Formulation  and  Model  Variants 

We  need  the  following  notation  (for  each  period  t)  to  describe  the  generic  single- 
item dynamic  lot  size  model. 
dt  = demand 
Xt  = production  level 

Dt  = cumulative  demand,  that  is,  Dt  = £U=1  dt 
Xt  = cumulative  production,  that  is,  Xt  = Xi 
Xt  = the  set  of  feasible  production  levels 
TZt  = the  set  of  feasible  cumulative  production  levels 
Pt(xt)  = cost  function  for  producing  xt  units 
Ht(Xt)  = cost  function  for  carrying  Xt  — Dt  units  inventory 
at  the  end  of  period  t 

Without  loss  of  generality,  we  assume  that  X0  = 0 ,D0  = 0 and  Xt  = Dj.  The 
generic  single-item  dynamic  lot  size  model  can  be  formulated  as  follows. 

(PI)  Z*  = min  + Ht(Xt)) 

(xu...,xT)  t=1 


such  that 


Xo  — 0, 

Xt  = Dt-, 

— Xt~  i + Xt, 
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xt  e xt,  v t 

Xt  € TZti  V t. 

Note  that  the  cost  functions,  Pt(xt ) and  as  well  as  the  feasible  sets,  Xt 

and  7 Zt,  in  the  above  model  are  not  specified.  By  specifying  these  functions  or  sets,  we 
can  define  variants  of  the  single-item  dynamic  lot  size  model.  We  need  the  following 
additional  notation  for  these  variants. 

Cf  = capacity  on  production  in  period  t. 
wt  = capacity  on  inventory  in  period  t. 

Ct  = cumulative  capacity  on  production,  that  is,  Ct  = J2l=i  ci- 
Kt  = fixed  setup  cost  in  period  t. 
pt  = unit  production  cost  in  period  t. 

= unit  purchasing  cost  in  period  t. 
ht  = unit  holding  cost  in  period  t. 
bt  = unit  backlogging  cost  in  period  t. 

S(x t)  = 0 if  xt  = 0,  or  1 if  xt  > 0. 

(Ml)  The  Wagner- Whitin  model 

In  the  Wagner- Whitin  model,  the  unit  production  cost  is  constant  ( pt  = p),  back- 
logging  is  not  allowed,  and  there  is  no  capacity  restraint  on  production  or  inventory. 
Thus  we  define 


Pt(xt)  = Kt8{xt)  + pxt 
Ht(Xt)  = ht(Xt-Dt) 
Xt  — [0,  oo) 

Pt  = [Dt,  Dt] 

(M2)  The  uncapacitated  dynamic  lot  size  model 
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The  uncapacitated  dynamic  lot  size  model  is  similar  to  the  Wagner- Whitin  model 
except  it  has  time  varying  unit  production  costs.  That  is, 

Pt(xt)  = KtS(xt ) + ptxt 

(M3)  The  dynamic  lot  size  model  with  capacity  on  production 

The  dynamic  lot  size  model  with  capacity  on  production  is  similar  to  (M2)  except 

Xt  = [0,  ct] 

Tit  = [Dt,mm{DT,Ct}] 

(M4)  The  dynamic  lot  size  model  with  capacity  on  inventory 

The  dynamic  lot  size  model  with  capacity  on  inventory  is  similar  to  (M2)  except 

Tit  = [Dt,mm{DT,wt  + Dt}] 

(Note  that  (M3)  and  (M4)  can  be  combined  to  have  capacities  on  both  produc- 
tion and  inventory.) 

(M5)  The  dynamic  lot  size  model  with  backlogging 

In  models  with  backlogging,  the  costs  of  carrying  negative  inventories  are  in- 
cluded, and  the  cumulative  production  level,  Xt,  can  be  less  than  the  cumulative 
demand  in  the  previous  period.  For  example,  (M4)  with  backlogging  is  obtained  with 
the  following  definitions: 


Pt(xt)  = 

Ht(Xt ) = 
Xt  = 


Kt6(xt)  + ptxt 

f bt(Dt  - Xt)  if  Xt  < Dt 
\ ht(Xt  - Dt)  if  > Dt 

[0,  oo) 


7it  = [0,min{Dr,u;t  + Dt}] 


(M6)  The  multi-source  dynamic  lot  size  model 

The  multi-source  model  means  that  there  is  more  than  one  way  to  satisfy  de- 
mands. For  example,  we  may  produce  in-house  or  buy  from  outside  suppliers.  Or  we 
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may  have  several  production  lines  with  independent  setups  or  different  production 
costs. 

The  difference  between  the  models  with  a single  source  and  those  with  multiple 
sources  is  only  in  the  production  cost  function.  For  instance,  in  the  make-or-buy 
cases,  suppose  that  the  in-house  capacity  is  ct  in  period  i,  and  let  pu  f3t  be  the  unit 
production  cost  and  the  unit  purchasing  cost,  respectively.  There  is  a fixed  charge 
for  each  case:  say,  Kt  is  the  setup  cost  for  production  and  at  is  the  setup  cost  for 
purchasing  in  period  t.  Also,  flt  is  assumed  to  be  greater  than  p*,  and  there  exists  e*, 
a break-even  point  between  producing  and  purchasing,  that  is,  Kt  + ptet  = at  + /3tet 
and  et  < ct.  Then  the  production  cost  function  of  this  model  is 

( atS(xt)  + ptxt  if  0 < xt  < et 

Pt(xt)  = l I<t  + Ptxt  if  et  < xt  < ct 

{ Kt  + ptct  + ott  + /3t(xt  - ct)  if  ct  < xt. 

Notice  that  above  function  is  not  necessarily  convex  or  concave. 

(M7)  The  multi-warehouse  dynamic  lot  size  model  In  this  model,  the  inventory  is  held 
in  several  warehouses.  Each  warehouse  has  its  own  capacity,  and  a fixed  charge  must 
be  paid  if  it  is  used.  Now,  suppose  there  are  n identical  warehouses  with  capacity  wt 
and  setup  cost  Wt  in  period  t.  Then  the  multi-warehouse  dynamic  lot  size  model  has 
the  following  holding  cost  function. 


Ht(Xt)  = 


0 if  Xt  - Dt  = 0 

iWt  + ht(Xt  — Dt)  if  ( i — 1 )wt  < Xt  — Dt  < iwt , 

i = 1 , . . . , n 


Notice  that  above  function  is  also  not  necessarily  convex  or  concave. 


3.1.2  A Continuous  Dynamic  Programming  Approach 

Let  Ft(Xt)  be  the  optimal  value  function  of  (PI)  considering  only  the  first  t 
periods  with  Xt  as  the  final  cumulative  production  level,  that  is, 

Ft(Xt)  = min  V (Pk(Xk)  + Hk{Xk)) 
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such  that 


*o  = 0 

Xk  X k — 1 T x ki  — 1 , • . . , 1 

% k G X A;  •>  k — • • • 5 ~t 

xkenk,  fc  = 1, . . . , t 1. 


And  define 


^o(*o)  = 

Then  (PI)  can  be  solved  by  dynamic 
tions: 


0 if  *0  = 0 

00  otherwise. 


programming  with  the  following  recursive  equa- 


Ft(Xt)  = Ht(Xt)  + min  {Ft_i  (**_!)  + Pt(xt) 


Xt  — Xt- 1 + Xt) 


Xt- 1 G Tlt-i)  xt  G V t 


Z*  = Fr(Ar)- 

The  conventional  approach  for  (3.1)  assumes  that  the  state  variable  A*  and  the 
control  variable  Xf  are  discrete.  Thus,  it  yields  a pseudo-polynomial  time  algorithm 
with  computational  complexity  O(DtX)  where  * is  the  number  of  total  elements 
in  Xt  over  all  t.  For  instance,  the  complexity  is  O(DtCt)  for  the  problems  in  (M3). 

Since  Dt  and  * can  be  very  large  numbers,  researchers  usually  use  this  approach 
only  as  a last  resort. 

Our  approach  is  different  from  the  conventional  one  by  treating  Xt  and  xt  as 
continuous  variables.  Since  Xt  is  continuous,  there  are  an  infinite  number  of  states 

in  our  approach.  In  order  to  resolve  the  difficulty  incurred  by  infinite  states,  the 
following  assumption  is  needed: 
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Assumption  Pt  (xt)  and  Ht(Xt)  are  piecewise  linear  functions  with  finitely  many 
segments  for  any  t. 

Since  all  given  cost  functions  are  piecewise  linear,  it  is  reasonable  to  expect  that 
I't, ( Xf ) is  also  a piecewise  linear  function  for  any  t.  The  following  theorem  will  be 
proven  later. 

Theorem  1 Ft(Xt)  is  a piecewise  linear  function  with  finitely  many  pieces  for  any 
t. 


Since  Ft-i(Xt-i)  is  piecewdse  linear,  its  domain  can  be  split  into  several  intervals, 


say 


ft.-.  = Uft!i. 

i 

such  that  in  each  interval  Ft-  is  only  a line  segment.  In  a similar  way, 

the  domain  of  Pfixf)  can  also  be  split  into  several  intervals,  that  is, 

Xt  = U xt(j\ 

3 

such  that  Pt(xt ) is  only  a single  line  segment  in  each  interval 

Now,  consider  minimizing  in  the  right  hand  side  of  (3.1)  segment  by  segment 
instead  of  point  by  point.  Define 


FtM(Xt)  = minjf^pQ-i)  + Pt(xt)  \ Xt  = Xt-x  + xu 


AVie^,  xt  e Xtu>}  V t,i,j. 


U) 


(3.2) 


Therefore,  the  recursive  equations  in  (3.1)  can  be  rewritten  as 


Ft(Xt)  = Ht(Xt)  + mm{Ft{i’j)(Xt)\\/i,j}  V t. 


The  function  value  of  F^(Xt)  is  infinite  if  Xt  is  out  of  its  domain. 


(3.3) 


In  the  following  sections  we  will  show  Ft  {^t)  f°r  any  ^ and  j is  a piecewise 
linear  function,  and  it  can  be  written  in  closed  form. 
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Given  a R-t-i  and  an  X}3\  suppose  that 

= \4\4\ 
x,{,)  = vP.tiX 

and 

Ft- i(*t-i)  = R^  + r^ Xt-u  for  Xt.x  G 

Pt(xt)  = P^+p^Xf  for  xt  G Xj3\ 

Then  for  each  i,  j,  (3.2)  becomes 

if^Xt)  = min  {f?w  + rW^_!  + Pb9  + p^x,  | 

^ = <4**^ -^t-i  < <4*  \ (3-4) 

We  use  a continuous  version  of  the  reaching  concept  (Denardo  (1982))  to  evaluate 
the  above  recursive  equation,  and  say  that  Xt  can  be  reached  from  Xt-i  if  b[3^  < 

Xt  - W-i  < b{j). 

Consider  the  minimization  in  (3.4).  Given  an  Xt,  after  substituting  Xt-\  = 
Xt  — Xt  or  xt  = Xt  — Xt-i,  we  are  to  minimize 

R{l)  - (rW  - p^)xt  + Pi3)  + r(l)Xt 
or 

R(l)  + (rw  - p(j))Xt- 1 + P{j)  + pU)Xt 

over  all  possible  production  levels  xt , or  over  all  Xt-\  from  which  Xt  can  be  reached. 
So,  if  rb)  < pb),  it  is  optimal  to  make  xt  as  small  as  possible,  or  to  make  Xt-\  as 
large  as  possible.  That  is,  xt  = b[3)  when  Xt  G [<4°  + b[j\a^  + b[j)]-  and  Xt-i  = af 
when  Xt  G [<4^  + b[3\  a^  + b^].  On  the  other  hand,  in  case  of  rb)  > pb)?  xt  should 
be  as  large  as  possible  and  Xt— i as  small  as  possible.  That  is,  Xt— i = when 
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V € [d\  ^ a{  ^ + &2^];  an(l  xt  — bP  when  Xt  € [<zj^  + <T  + 4^]-  The  above 

discussion  can  be  summarized  in  the  following  algebraic  expressions:  If  rb)  < p^\ 
then 


otherwise, 


< 


R(i)  - (r^  - p^)b[j)  + pV)  + r®Xt, 

for  aP  + bP  < Xt  < aip  + bP 


R(i)  + (r(d_pb))ab)  + p(i)  + p(i)xt, 

for  a ^ < Xt  < a,P  + b[3\ 


/ 


< 


V 


/?(*)  + (rb')-pb))a[!)  + pO)  + pO')^, 

for  _).  bti)  < Xt  < af*  + bP 

R(i)  _ (rb)_p0))^)  + p(j)  + r^Xt, 

for  ap  + hip  < Xt  < aP  + bp, 


or  in  another  form,  if  rb)  < pb)?  then 


otherwise, 


p{%j) 


(*<) 


' F<- i(x,  - ti’1)  + p,(^\ 

for  o,!'  + h\Jt  < X,  < 

Tt-i(<4  + Pt{Xt  — <4  ^)> 

k for  <4*'  + &P  < < aip  + 6^, 


' Ft^(aP)  + Pt(Xt-a{p), 

for  aP  + bP  < Xt  < ap  + bP 

F,-i(X,  - bf)  + 

, for  + /ojl  < Xt<  a,’  + /o' 1 . 


We  interpret  the  above  segments  as  follows.  The  first  segment  of  (3.5)  and  the 
second  segment  of  (3.6)  are  created  by  fixing  the  production  levels  at  b[J^  and  b[J\ 
respectively.  On  the  other  hand,  the  second  segment  of  (3.5)  and  the  first  segment  of 
(3.6)  are  created  by  fixing  the  cumulative  production  level  in  period  t — 1,  Xt- 1,  at 
and  a±  \ respectively. 

We  have  shown  that  F}  ^\Xt)  in  (3.3)  for  any  i and  j is  a piecewise  linear 
function  with  two  segments.  Therefore  Ft(Xt)  in  (3.3)  can  be  computed  in  two  steps: 
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(c)  (d) 

Figure  3.1.  Segment  update  of 

first,  evaluating  the  minimum  of  over  all  i and  j,  then  adding  this  result  to 

Ht(Xt). 

The  hrst  step  is  illustrated  in  Figure  3.1.  Let  there  be  a total  of  s line  segments 
from  all  of  the  F1/  ’^\Xt)  over  all  i and  j.  Denote  L ^ as  the  &th  line  segment  for 
1 < k < s (taken  in  any  order).  Any  L\  is  uniquely  defined  by  ( l,u,p,y ) = (lower 
bound,  upper  bound,  slope,  intercept);  see  Figure  3.1(b). 

Starting  with  Ft(1)(Xt)  = L\  for  Xt  € [l,u]f)Kt,  the  function  Ft(s)(Xt)  for 
Xt  € F-t  can  be  obtained  by  a segment  by  segment  update  procedure.  In  Fig- 
ure 3.1(a),  FF~1]  represents  Ft  after  (k  - 1)  updates.  It  is  defined  by  its  breakpoints 
X(0),  X(l),  A'(2),X(3)  and  by  its  segments  which  have  slopes  P(l),  P(2),  P(3)  and 
intercepts  F(l),  Y(2),  F(3).  Ft{k~1]  and  L\  are  both  included  in  Figure  3.1(c),  then 
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(k)  • 

Ft  is  obtained  by  moving  through  the  breakpoints,  including  / and  u,  and  comparing 
the  function  values  with  those  of  the  segment.  The  result  is  given  in  Figure  3.1(d). 

The  second  step  is  to  sum  two  piecewise  linear  functions,  and  F}s\xt). 

After  sorting  all  the  breakpoints  of  these  two  functions,  there  is  only  a single  line 
segment  for  each  function  between  two  such  consecutive  breakpoints.  Therefore,  the 
summation  of  these  two  functions  can  be  calculated  easily  to  give  Ft(Xt). 

For  a detailed  discussion  ol  evaluating  lower  envelopes  and  recovering  an  optimal 
solution  for  the  model  (M3),  see  Chen  et  al.  (1992),  which  contains  a pseudocode 
and  a numerical  example. 

The  computational  effort  of  our  approach  is  a function  of  the  time  horizon, 
T,  the  number  of  segments  in  the  cost  functions  and  the  number  of  segments  in 
the  optimal  value  functions.  Although  theoretically  the  number  of  segments  may 
increase  at  an  exponential  rate  with  T,  it  need  not  be  more  than  Dj  in  practice  since 
we  can  always  delete  segments  between  two  consecutive  integer  values  of  Xt  without 
losing  the  optimal  solution.  So,  our  approach  can  not  be  worse  than  the  conventional 
discrete  approach  for  problems  with  integer  data. 

We  conclude  this  section  by  proving  Theorem  1. 

Proof  of  Theorem  1 The  theorem  will  be  proven  by  mathematical  induction. 
When  t = 0,  F0(Xo)  is  a piecewise  linear  function  with  one  segment  (the  original 
point)  by  definition.  Now  assume  that  Ft_i(Xt-i)  is  a piecewise  linear  function  with 
finite  segments  for  t > 1.  From  (3.3),  formulas  (3.5),  (3.6)  and  the  assumption  that 
Ht(Xt)  is  a piecewise  linear  function  with  finite  segments,  it  can  be  seen  that  Ft(Xt) 
is  also  a piecewise  linear  function  with  finite  segments.  | 

3.1.3  Computational  Experience 

In  this  section  we  examine  the  efficiency  of  our  approach  on  the  capacitated 
DLSM  with  multiple  piece  production  cost  functions.  Suppose  there  are  M pieces 
in  the  production  cost  function  of  each  period.  There  is  a fixed  charge  (setup  cost), 
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a slope  (unit  production  cost),  and  a capacity  (length  of  such  piece)  for  each  piece. 
Linear  holding  costs  are  considered  in  this  model  and  no  backlogging  is  allowed. 

We  modified  the  problem  pattern  in  Baker  et  al.  (1978)  to  randomly  create  test 
problems. 

The  demand  is  given  in  a pattern  as  follows: 

2 7 r 

dt  = 200  + azt  + a sin  [— (t  + b/4)]  (3.7) 

where  a = standard  error  of  demand, 

zt  = i.i.d.  standard  normal  random  deviates, 
a = amplitude  of  the  seasonality  component, 
b = length  of  seasonal  cycle  in  periods. 

Five  problems  were  created  in  each  of  the  following  four  combinations:  (1)  a = 

67,  a = 0,  (2)  (j  = 237,  a = 0,  (3)  a = 67,  a = 125,  b = T and  (4)  a = 67,  a = 125, 
6=  12. 

We  tested  the  problems  with  T € {24, 48, 96, 192}  and  M G {1, 2, 4, 8}.  The  unit 
holding  cost,  ht,  for  each  t is  uniformly  distributed  between  0.5  and  1.5.  The  unit  pro- 
duction cost,  Pt  , for  each  period  t and  each  piece  j is  uniformly  distributed  between 
10  and  30.  The  capacity  for  each  period  and  each  piece  is  uniformly  distributed  in 
(0.5C/M,  1.5C/M)  where  C 6 {400,800,1200,1600}.  These  average  total  capacities 
C satisfy  (average)  demand  of  2,  4,  6 and  8 periods,  respectively.  The  correspond- 
ing setup  costs  are  K = 400, 1600,3600,6400  which  come  from  the  continuous  EOQ 
(economic  order  quantity)  model.  The  setup  cost  for  each  period  and  each  piece  is 
uniformly  distributed  in  (0.5 K/M,  1.5 K/M). 

A total  of  5120  feasible  problems  were  created  and  tested.  The  program  of  our 
algorithm  is  coded  in  FORTRAN  and  run  on  a VAX  8600. 

We  are  interested  in  both  cpu  time  and  the  number  of  segments  of  the  optimal 
value  function.  The  results  are  listed  in  two  tables:  one  for  cpu  time  and  the  other 
for  the  maximal  number  of  segments.  Since  these  tables  are  large,  they  are  placed  at 
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the  end  of  this  section  (Table  3.5  to  Table  3.8).  Each  cell  in  these  tables  represents 
20  different  demands  of  the  corresponding  category.  Both  average  and  the  maximal 
(in  the  parentheses)  values  of  the  20  problems  are  listed.  Note  that  the  tables  show 

that  these  computational  results  are  very  stable  since  the  maximal  values  are  within 
twice  the  average  values  in  almost  every  case. 

As  expected,  the  number  of  segments  of  the  optimal  value  function  is  very  small. 
Even  though  this  number  can  be  more  than  2T  for  a T-period  problem,  it  is  only  a 
fraction  of  its  theoretical  value  as  shown  in  Table  3.1.  The  fact  that  the  number  of 
segments  is  small  ensures  that  problems  as  large  as  96  periods  and  8 segments  in  the 
production  cost  function  can  be  solved  in  less  than  100K  bytes  of  memory  since  about 
10  bytes  are  required  for  each. 


Table  3.1  Maximal  segments  in  Ft(Xt)  over  time  horizon 


T 

12 

24 

48 

96 

21 

4.1  x 103 

1.7  x 107 

2.8  x 1014 

7.9  x 1028 

actual  number 

607 

1704 

2988 

6102 

The  results  in  the  Appendix  also  show  that  the  problems  with  smaller  setup 
costs  are  easier  for  our  algorithm  than  those  with  larger  setup  costs.  On  the  other 
hand,  the  problems  with  less  capacity  are  harder  than  those  with  more  capacity.  Out 
of  the  16  possible  combinations  of  K and  C,  the  case  with  K = 6400  and  C = 400  is 
the  hardest.  Table  3.2  lists  the  average  cpu  time  of  this  hardest  case.  In  the  case  of 
T = 96  and  M = 1,  it  takes  only  1.97  seconds  on  average  to  find  the  optimal  solution. 
We  feel  these  results  demonstrate  the  effectiveness  of  the  method  relative  to  the  prior 
literature.  For  example,  the  largest  problem  previously  solved  to  optimality,  with 
M = 1,  only  has  T = 18  (Kirca  (1990)). 


Table  3.2  Average  cpu  time  (K=6400,  C=400) 


T =12 

T = 24 

oo 

T = 96 

M = 1 

0.01 

0.06 

0.40 

1.97 

M — 2 

0.05 

0.27 

1.84 

12.49 

M = 4 

0.19 

1.59 

9.95 

67.76 

M = 8 

1.38 

10.89 

75.69 

512.44 

Note:  seconds  on 

Vax  8600 
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Finally,  the  increases  in  the  cpu  times  in  Table  3.2  are  given  in  Table  3.3  and 
Table  3.4  as  a function  of  T and  M , respectively. 


Table  3.3  Increases  in  cpu  time  with  respect  to  T 


T 

12  -»•  24 

24  48 

48  ->  96 

M = 1 

6.0 

6.7 

4.9 

M = 2 

5.4 

6.8 

6.8 

M = 4 

8.4 

6.3 

6.8 

M = 8 

7.9 

7.0 

6.8 

Table  3.4  Increases  in  cpu  time  with  respect  to  M 


M 

T=  12 

T = 24 

II 

oo 

to 

05 

II 

1 2 

5.0 

4.5 

4.6 

6.3 

2 — ► 4 

3.8 

5.9 

5.4 

5.4 

4 — ► 8 

7.3 

6.8 

7.6 

7.6 

In  these  tables  an  entry  of  8 indicates  a cubic  growth  rate.  The  entries  in  the 
two  tables  range  from  3.8  to  8.4,  demonstrating  that  the  growth  rates  range  from  1.93 
to  3.07.  Thus  it  is  reasonable  to  claim  that,  for  data  similar  to  ours,  the  cpu  time 
increases  at  a rate  which  is  approximately  cubic  in  either  T or  M. 
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Table  3.5  Cpu  time  (seconds  on  a VAX  8600) 


T = 12 

T = 24 

K 

C 

M=  1 

M = 2 

M = 4 

M = 8 

M=  1 

M=  2 

M=  4 

M=  8 

400 

400 

0.01“ 

0.02 

0.09 

0.33 

0.03 

0.11 

0.39 

1.68 

(0.01)6 

(0.04) 

(0.12) 

(0.49) 

(0.05) 

(0.14) 

(0.54) 

(2.37) 

800 

0.01 

0.01 

0.05 

0.17 

0.02 

0.06 

0.20 

0.83 

(0.01) 

(0.03) 

(0.07) 

(0.23) 

(0.03) 

(0.07) 

(0.24) 

(1.05) 

1200 

0.00 

0.01 

0.03 

0.11 

0.01 

0.04 

0.14 

0.51 

(0.01) 

(0.02) 

(0.05) 

(0.16) 

(0.02) 

(0.06) 

(0.19) 

(0.67) 

1600 

0.00 

0.01 

0.02 

0.08 

0.01 

0.03 

0.10 

0.36  i 

(0.01) 

(0.02) 

(0.04) 

(0.10) 

(0.01) 

(0.04) 

(0.14) 

(0.46) 

1600 

400 

0.01 

0.03 

0.13 

0.68 

0.04 

0.18 

0.83 

4.66 

(0.02) 

(0.05) 

(0.25) 

(1.16) 

(0.06) 

(0.26) 

(1.19) 

(7.03) 

800 

0.00 

0.02 

0.07 

0.32 

0.02 

0.09 

0.42 

1.96 

(0.01) 

(0.04) 

(0.11) 

(0.48) 

(0.04) 

(0.13) 

(0.67) 

(2.52) 

1200 

0.00 

0.01 

0.05 

0.18 

0.01 

0.05 

0.22 

1.04 

(0.02) 

(0.02) 

(0.06) 

(0.27) 

(0.02) 

(0.08) 

(0.31) 

(1.66) 

1600 

0.01 

0.01 

0.03 

0.12 

0.01 

0.04 

0.13 

0.62 

(0.01) 

(0.02) 

(0.06) 

(0.17) 

(0.02) 

(0.06) 

(0.18) 

(0.81) 

3600 

400 

0.01 

0.04 

0.20 

0.90 

0.05 

0.23 

1.03 

8.39 

(0.02) 

(0.06) 

(0.33) 

(1.51) 

(0.08) 

(0.36) 

(1.84) 

(18.87) 

800 

0.01 

0.02 

0.09 

0.37 

0.03 

0.13 

0.60 

3.64 

(0.02) 

(0.04) 

(0.18) 

(0.60) 

(0.05) 

(0.20) 

(0.93) 

(6.37) 

1200 

0.01 

0.02 

0.06 

0.22 

0.02 

0.07 

0.35 

1.46 

(0.02) 

(0.03) 

(0.10) 

(0.30) 

(0.03) 

(0.11) 

(0.49) 

(2.00) 

1600 

0.00 

0.01 

0.04 

0.15 

0.01 

0.05 

0.20 

0.90 

(0.01) 

(0.02) 

(0.07) 

(0.24) 

(0.03) 

(0.08) 

(0.29) 

(1.29) 

6400 

400 

0.01 

0.05 

0.19 

1.38 

0.06 

0.27 

1.59 

10.89 

(0.02) 

(0.09) 

(0.34) 

(3.93) 

(0.08) 

(0.44) 

(2.71) 

(39.42) 

800 

0.01 

0.03 

0.10 

0.54 

0.04 

0.16 

0.80 

4.10 

(0.02) 

(0.04) 

(0.16) 

(l.n) 

(0.06) 

(0.27) 

(1.41) 

(5.97) 

1200 

0.01 

0.01 

0.06 

0.29 

0.02 

0.09 

0.42 

2.06 

(0.02) 

(0.03) 

(0.10) 

(0.43) 

(0.04) 

(0.17) 

(0.70) 

(3.62) 

1600 

0.00 

0.01 

0.04 

0.16 

0.02 

0.06 

0.28 

1.29 

(0.01) 

(0.02) 

(0.07) 

(0.25) 

(0.03) 

(0.08) 

(0.43) 

(1.98) 

a : average  b : worst  case 
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Table  3.6  Cpu  time  (seconds  on  a VAX  8600) 


T = 12 

T = 24 

K 

C 

M—  1 

M=  2 

M=  4 

M=  8 

1 

M=  2 

M=  4 

M—  8 

400 

400 

0.11“ 

0.41 

1.64 

7.15 

0.42 

1.64 

7.26 

32.72 

(0.13)6 

(0.51) 

(1.94) 

(8.94) 

(0.55) 

(1.99) 

(8.11) 

(38.85) 

800 

0.07 

0.25 

0.89 

3.73 

0.29 

1.05 

3.97 

16.47 

(0.10) 

(0.35) 

(1.07) 

(4.29) 

(0.33) 

(1.22) 

(4.86) 

(20.30) 

1200 

0.04 

0.16 

0.59 

2.23 

0.18 

0.68 

2.59 

9.92 

(0.09) 

(0.21) 

(0.78) 

(2.85) 

(0.22) 

(0.90) 

(3.05) 

(12.87) 

1600 

0.04 

0.12 

0.43 

1.67 

0.13 

0.47 

1.80 

6.95 

(0.07) 

(0.16) 

(0.52) 

(1.96) 

(0.18) 

(0.55) 

(2.56) 

(7.93) 

1600 

400 

0.19 

0.82 

4.38 

28.27 

0.79 

3.60 

20.15 

128.24 

(0.31) 

(1.01) 

(6.16) 

(44.79) 

(1.07) 

(4.92) 

(26.19) 

(171.49) 

800 

0.11 

0.41 

1.73 

9.47 

0.44 

1.72 

8.31 

46.71 

(0.16) 

(0.55) 

(2.17) 

(13.70) 

(0.51) 

(2.28) 

(10.23) 

(63.35) 

1200 

0.06 

0.24 

0.99 

4.54 

0.26 

1.04 

4.22 

22.84 

(0.09) 

(0.29) 

(1.17) 

(5.88) 

(0.34) 

(1.23) 

(5.01) 

(33.48) 

1600 

0.05 

0.17 

0.67 

2.95 

0.19 

0.70 

2.94 

13.93 

(0.08) 

(0.23) 

(0.80) 

(3.92) 

(0.25) 

(0.84) 

(3.51) 

(18.09) 

3600 

400 

0.28 

1.24 

8.38 

51.62 

1.37 

7.18 

40.48 

334.90 

(0.41) 

(1.70) 

(13.07) 

(102.52) 

(1.89) 

(9.51) 

(58.75) 

(559.83) 

800 

0.15 

0.65 

3.52 

19.91 

0.65 

3.04 

16.96 

119.41 

(0.23) 

(0.90) 

(5.76) 

(29.38) 

(0.78) 

(3.63) 

(26.28) 

(156.69) 

1200 

0.08 

0.36 

1.61 

10.90 

0.37 

1.64 

7.74 

51.35 

(0.13) 

(0.48) 

(1.96) 

(14.73) 

(0.47) 

(2.37) 

(9.33) 

(63.94) 

1600 

0.05 

0.23 

0.99 

5.35 

0.25 

1.02 

4.67 

28.29 

(0.07) 

(0.34) 

(1.26) 

(7.80) 

(0.31) 

(1.24) 

(5.54) 

(36.93) 

6400 

400 

0.40 

1.84 

9.95 

75.69 

1.97 

12.49 

67.76 

512.44 

(0.65) 

(2.54) 

(23.88) 

(145.98) 

(2.91) 

(27.91) 

(94.17) 

(766.11) 

800 

0.19 

0.82 

4.98 

39.08 

0.93 

4.88 

31.11 

210.95 

(0.27) 

(1.07) 

(7.14) 

(59.82) 

(1.19) 

(7.73) 

(45.10) 

(281.85) 

1200 

0.10 

0.48 

2.25 

18.73 

0.50 

2.42 

12.43 

88.28 

(0.16) 

(0.66) 

(2.84) 

(38.88) 

(0.69) 

(3.14) 

(17.62) 

(112.00) 

1600 

0.07 

0.30 

1.48 

8.68 

0.30 

1.45 

7.76 

50.86 

(0.10) 

(0.47) 

(2.34) 

(11.49) 

(0.40) 

(1.87) 

(9.50) 

(62.02) 

a : average  b : worst  case 
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Table  3.7  Maximal  segments  of  Ft(Xt)  over  time  horizon 


K 

C 

II 

oo 

1 

T = 96 

M=  1 

M=  2 

II 

M=  8 

M=  1 

CM 

II 

II 

M — 8 

400 

400 

12.4“ 

26.9 

55.1 

102.2 

22.6 

53.9 

106.4 

216.3 

(22)6 

(41) 

(79) 

(133) 

(34) 

(74) 

(144) 

(271) 

800 

8.4 

16.1 

31.0 

56.5 

15.9 

30.5 

59.3 

118.9 

(13) 

(21) 

(41) 

(85) 

(23) 

(43) 

(74) 

(144) 

1200 

6.4 

11.9 

21.7 

38.5 

11.1 

22.0 

43.8 

79.8 

(9) 

(17) 

(29) 

(55) 

(14) 

(27) 

(66) 

(103) 

1600 

5.6 

8.8 

15.6 

28.3 

10.1 

17.1 

32.2 

58.5 

(8) 

(12) 

(22) 

(35) 

(14) 

(21) 

(42) 

(77) 

1600 

400 

15.8 

38.1 

81.3 

178.8 

33.5 

83.2 

203.9 

429.8 

(21) 

(63) 

(124) 

(273) 

(59) 

(121) 

(295) 

(565) 

800 

8.8 

19.0 

40.8 

94.7 

21.3 

45.0 

115.6 

217.8 

(ii) 

(26) 

(62) 

(131) 

(35) 

(74) 

(176) 

(300) 

1200 

7.7 

13.8 

30.0 

56.7 

9.6 

30.1 

65.1 

136.7 

(ii) 

(19) 

(44) 

(86) 

(20) 

(46) 

(92) 

(220) 

1600 

5.9 

10.6 

21.1 

37.2 

11.3 

23.0 

45.7 

92.0 

(8) 

(15) 

(40) 

(53) 

(17) 

(34) 

(73) 

(117) 

3600 

400 

16.7 

41.0 

114.8 

217.0 

43.7 

110.4 

244.9 

608.7 

(25) 

(74) 

(185) 

(308) 

(76) 

(191) 

(401) 

(1000) 

800 

10.4 

22.7 

56.0 

101.3 

26.5 

61.0 

151.6 

342.6 

(15) 

(36) 

(107) 

(153) 

(42) 

(133) 

(191) 

(575) 

1200 

7.5 

16.6 

34.3 

69.0 

16.6 

37.0 

98.6 

178.1 

(14) 

(28) 

(60) 

(95) 

(27) 

(54) 

(141) 

(251) 

1600 

5.9 

11.6 

25.1 

49.5 

12.4 

28.8 

60.3 

122.9 

(9) 

(24) 

(53) 

(71) 

(18) 

(42) 

(79) 

(183) 

6400 

400 

17.9 

48.3 

107.5 

293.6 

49.8 

129.3 

337.1 

714.2 

(28) 

(79) 

(177) 

(607) 

(83) 

(207) 

(580) 

(1704) 

800 

10.4 

24.0 

58.3 

134.7 

28.9 

73.3 

178.5 

350.3 

(16) 

(40) 

(90) 

(229) 

(54) 

(109) 

(288) 

(436) 

1200 

7.5 

16.4 

35.7 

90.1 

18.3 

42.7 

108.5 

222.1 

(12) 

(25) 

(56) 

(186) 

(28) 

(81) 

(167) 

(343) 

1600 

5.7 

11.9 

26.3 

53.8 

14.1 

27.6 

77.0 

159.3 

(8) 

(18) 

(45) 

(86) 

(21) 

(36) 

(127) 

(256) 

a : average  b : worst  case 
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Table  3.8  Maximal  segments  of  Ft(Xt)  over  time  horizon 


l 

00 

II 

bs 

T = 96 

K 

C 

M=  1 

M—  2 

II 

*5 

M=8 

M—  1 

M = 2 

II 

M = 8 

400 

400 

40.5“ 

98.2 

212.1 

405.6 

83.5 

190.8 

428.8 

841.4 

(57)‘ 

(129) 

(253) 

(513) 

(117) 

(234) 

(491) 

(937) 

800 

28.1 

64.2 

118.9 

235.1 

57.0 

118.9 

244.8 

462.0 

(35) 

(107) 

(156) 

(295) 

(74) 

(142) 

(275) 

(571) 

1200 

20.4 

43.0 

83.9 

158.3 

37.0 

81.8 

168.5 

303.8 

(31) 

(57) 

(in) 

(196) 

(47) 

(99) 

(193) 

(376) 

1600 

16.7 

32.5 

63.7 

128.0 

28.8 

60.7 

124.7 

228.1 

(23) 

(50) 

(83) 

(163) 

(36) 

(73) 

(162) 

(264) 

1600 

400 

78.5 

191.2 

460.6 

1037.7 

153.3 

389.5 

937.7 

2115.8 

(124) 

(250) 

(635) 

(1562) 

(209) 

(508) 

(1167) 

(2846) 

800 

41.0 

102.7 

216.1 

449.0 

85.7 

199.6 

439.8 

936.0 

(56) 

(152) 

(273) 

(595) 

(in) 

(269) 

(641) 

(1173) 

1200 

28.5 

61.5 

127.3 

280.5 

58.0 

121.8 

265.0 

559.0 

(48) 

(77) 

(169) 

(348) 

(76) 

(151) 

(322) 

(764) 

1600 

21.3 

47.0 

97.1 

203.1 

44.3 

87.2 

196.9 

372.7 

(38) 

(70) 

(113) 

(276) 

(61) 

(103) 

(244) 

(462) 

3600 

400 

102.7 

271.8 

727.4 

1494.4 

243.3 

682.0 

1537.9 

3692.5 

(150) 

(409) 

(941) 

(2323) 

(356) 

(916) 

(2179) 

(4879) 

800 

56.0 

140.8 

374.5 

733.8 

121.3 

329.9 

775.2 

1793.6 

(90) 

(191) 

(553) 

(988) 

(151) 

(419) 

(1129) 

(2157) 

1200 

35.5 

91.3 

206.3 

490.2 

75.7 

193.4 

430.6 

947.3 

(57) 

(128) 

(248) 

(633) 

(101) 

(270) 

(528) 

(1208) 

1600 

25.3 

61.3 

139.4 

310.6 

57.2 

124.7 

289.7 

639.6 

(43) 

(84) 

(200) 

(419) 

(70) 

(166) 

(372) 

(831) 

6400 

400 

156.6 

384.0 

795.0 

1922.3 

363.0 

1046.9 

2230.0 

4483.7 

(262) 

(530) 

(1435) 

(2988) 

(516) 

(2039) 

(2944) 

(6102) 

800 

75.2 

176.9 

463.1 

1113.1 

179.4 

474.0 

1145.6 

2533.0 

(105) 

(239) 

(616) 

(1586) 

(224) 

(691) 

(1458) 

(3671) 

1200 

44.8 

112.6 

257.9 

641.0 

101.3 

260.0 

582.4 

1401.8 

(68) 

(139) 

(337) 

(984) 

(167) 

(359) 

(830) 

(1744) 

1600 

33.5 

74.7 

193.1 

389.1 

68.4 

180.4 

411.1 

968.1 

(54) 

(103) 

(258) 

(531) 

(90) 

(232) 

(520) 

(1251) 

a : average  b : worst  case 


43 


3.2  Models  with  Two  Control  Variables 

In  the  first  three  sections  of  this  chapter,  the  models  we  discuss  have  only  one 
control  variable,  the  production  level  xt , as  well  as  one  corresponding  state  variable, 
the  cumulative  production  level  Xt . In  this  section,  we  extend  the  algorithm  in  Sec- 
tion 3.1  to  solve  models  with  two  control  variables  and  two  state  variables.  The  new 
state  variable  is  related  to  machine  status,  for  example,  machine  age,  total  number 
of  setups,  or  whether  the  machine  is  on  or  not.  The  new  control  variable  changes 
the  status  of  the  machine.  We  only  consider  cases  where  the  new  state  variable  is 
discrete. 

3.2.1  Formulation 

We  need  the  following  notation  to  describe  the  generic  single-item  dynamic  lot 
size  model  with  two  control  variables. 

dt  = demand  in  period  t. 

xt  — production  level  in  period  t. 

Dt  = cumulative  demand,  that  is,  Dt  = d{. 

Xt  = cumulative  production,  that  is,  Xt  = x% ■. 

xjt  = control  variable  for  machine  status  in  period  t. 

Yt  = machine  status  in  period  t. 

Xt(Yt-\,yt)  = a set  of  feasible  production  levels,  that  is,  xt  £ 

dlt  = a set  of  feasible  cumulative  production  levels,  that  is,  Xt  £ Kt. 

Yt  = a set  of  feasible  yt,  that  is,  yt  £ yt. 

St  = a set  of  feasible  Yt,  that  is,Ti  £ St. 

Pt{xtiYt-i,yt)  = cost  function  for  producing  xt  units  in  period  t. 

Ht(Xt)  = cost  of  carrying  Xt  — Dt  units  inventory  at  the  end  of  period  t. 

Note  that  Xt{Yt-i,yt)  as  well  as  Pt(xt,Yt-i,yt)  are  functions  of  the  machine 
status  Yt- 1 and  the  action  to  change  this  status  yt.  This  will  properly  reflect  the  real 
situation.  In  the  real  world,  the  production  rate  (and  thus  the  production  capacity) 
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c4  a machine  depends  on  the  machine’s  age  or  other  status.  For  instance,  an  old 
machine  has  a greater  possibility  to  break  down  then  a new  one.  Therefore,  an  old 
machine  has  less  capacity  than  a new  one. 

Without  loss  of  generality,  we  may  assume  that  X0  — 0,  Y0  = 0,  and  Xj  — Dj, 
and  for  notational  convenience,  let  Do  = 0.  The  generic  DLSM  with  two  control 
variables  can  be  formulated  as  follows. 

(P2)  Z*  = min  Yt_uyt)  + Ht{Xt)) 

* = 1 

such  that 

Xt  = Xt~ i + x-t , V t 
Yt  = g(Yt-Uyt),  Vi 

xt  € Xt{Yt~i,yt),  Vi 

yt  eyt,  Vi 

A t 6 P-ti  V i 

Yt  € St,  Vi. 

where  g(Yt-i,yt)  is  a transition  function  of  the  state  variable  Yt. 

Note  that  the  cost  functions,  Pt(xt^Yt-i,yt)  and  Ht{Xt ),  as  well  as  the  feasible 
sets,  Xt{Yt-x,yt),  Xt,  St , and  the  transition  function  g(Yt-i,yt)  in  the  above  model 
are  not  specified.  By  specifying  these  functions  or  sets,  we  may  define  variants  of  the 
single-item  dynamic  lot  size  model.  We  need  the  following  additional  notation  for  the 
description  of  these  variants. 

ct  = capacity  on  production  in  period  t. 

Ct  = cumulative  capacity  on  production,  that  is,  Ct  = Y?i=i  ci- 
Kt  — fixed  setup  cost  in  period  t. 
pt  = unit  production  cost  in  period  t. 
ht  = unit  holding  cost  in  period  t. 
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(M7)  The  DLSM  with  startup  and  reservation  costs 

The  DLSM  with  startup  and  reservation  costs  allows  two  kinds  of  fixed  costs. 
A startup  cost  is  charged  whenever  the  machine  is  switched  from  “off”  to  “on”.  A 
reservation  cost  is  incurred  in  any  period  when  the  machine  is  on,  whether  there  is 
production  or  not. 

There  are  two  states  for  the  second  state  variable  Yt  in  this  model:  Yt  — 0 means 
that  machine  is  “off”  and  Yt  = 1 means  it  is  “on” . There  are  also  two  possible  actions 
to  change  the  machine  status:  yt  — 0 means  “turn-machine-off”  for  period  t and 
yt  = 1 means  “turn-machine-on”  for  period  t. 

Let  rt  and  St  be  the  reservation  cost  and  startup  cost  in  period  t,  respectively, 
then  the  formulation  of  the  DLSM  with  reservation  and  startup  costs  is  (P2)  with 
the  following  definitions  of  functions  and  sets. 


Pt(xtiYt—hyt)  — 

Ht(Xt)  = 
g(Yt-i,yt)  = 

Xt(Yt-i,yt ) = 

yt  = 


{ o if  yt  = o 

< rt  + Ptxt  if  yt  - 1 and  V)_i  = 1 
{ st  + ptxt  if  yt  = 1 and  Yt- 1 = 0 

ht(Xt  - Dt) 

f 0 if  yt  = 0 
\ 1 if  Vt  = 1 

f {0}  if  yt  = 0 
1 (0,  ct]  if  yt  = 1 

{0,1} 


Ht  = [Dt,nun{Ct,DT}] 

St  = {0,1} 


(M8)  The  DLSM  with  learning  in  setups 

A DLSM  with  learning  in  setups  is  DLSM  in  which  setup  cost  is  a function  of 
the  total  number  of  setups.  This  model  comes  from  the  fact  that  workers  can  gain 
experience  (learn)  from  previous  set  ups,  thus  the  setup  time  as  well  as  setup  costs 
will  decrease  as  the  number  of  setups  increases. 
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In  this  model,  the  second  state  variable  Yt  represents  the  total  number  of  set 
ups  from  period  1 to  period  t.  There  are  two  possible  values  of  yt.  0 or  1.  The  former 
means  not  to  set  up  the  machine,  and  the  later  means  to  set  up. 

Let  Kn  be  the  setup  cost  of  the  nth  set  up,  then  the  formulation  of  the  DLSM 
with  learning  in  setups  is  (P2)  with  the  following  definitions  of  functions  and  sets. 


Yt— i , yt) 

Ht(Xt) 

g(Yt-i,yt) 

Xt(Y*-uyt) 

y* 

nt 

St 


I 0 if  yt  = 0 

\ Kn+ 1 + ptxt  if  yt  = 1 and  Yt- 1 = n 

= ht(Xt-Dt) 

= ( Yt- 1 if  yt  = 0 
\ Yt- 1 + 1 if  yt  = 1 

= f {0}  if  yt  = o 
\ (0,  ct]  if  yt  = 1 

= {0,1} 

= [Dt,mm{Ct,DT}] 

= {0 


(M9)  The  DLSM  with  machine  replacement 

In  the  DLSM  with  machine  replacement,  the  age  of  a machine  increases  when  it 
is  set  up  in  a production  period,  and  the  production  cost  as  well  as  the  capacity  are 
functions  of  the  machine’s  age.  In  each  period  we  are  to  decide  whether  to  replace  a 
new  machine  or  how  many  units  to  produce.  It  is  assumed  that  the  machine  is  always 
set  up  in  the  period  in  which  this  machine  is  replaced. 

In  this  model  the  second  state  variable  Yt  is  the  machine’s  age.  The  correspond- 
ing control  variable  yt  is  not  to  produce  anything  ( yt  = 0),  to  set  up  a machine 
(yt  = 1),  or  to  replace  a machine  (y\  = 2). 

Let  AT,  pn , cn,  and  rn  be  the  setup  cost,  unit  production  cost,  production 
capacity,  and  machine  residual  value  for  the  machine  at  age  n.  Let  Qt  be  the  price 
of  a new  machine  at  period  t.  The  formulation  of  the  DLSM  with  learning  in  setups 
is  (P2)  with  the  following  definitions  of  functions  and  sets. 
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Pt(xt,Yt-i,yt)  = 
Ht{Xt)  = 


g(Yt-i,yt)  = 


0 

< Fn+1  + Pn+lxt 

, Qt~rn  + K\  + pixt 
ht{Xt  - Dt) 


< 


+ l 

1 


if  yt  = o 
if  y t = l 
if  yt  = 2 


if  yt  = 0 

if  yt  = 1 and  Ft_i  = n 
if  yt  — 2 and  >7-i  — n 


( {0}  if  yt  = 0 

< (0,  cn+1]  if  yt  = 1 and  Yt-\  = n 

l (0,  ci]  if  = 2 

{0,1,2} 


77-t  = [Z>i,min{f  x ci,Dt}] 
= {0, 


3.2.2  Algorithm 


Similar  to  Section  3.1,  it  can  be  seen  that  (P2)  can  be  solved  by  our  dynamic 
programming  algorithm.  Let  Ft(Xt,Yt)  be  the  optimal  value  function  of  (P2)  consid- 
ering only  the  first  t periods  with  Xt  as  the  final  cumulative  production  level  and  Yt 
as  the  machine  status,  that  is, 


t 

Ft(Xt,  Yt)  = min  ^ (Pk(xk,  Yt-i,yk ) + Hk(Xk)) 

k=l 

such  that 


Xk  ^k—1  + %ki  k — I5  • • • 5 t 

Yk  = g{Yk-uyk),  k = l,...,t 

G %k{Yt— 1 iDk)-)  k = 1, . . . , t 


Uk  G 34,  k = 1, . . . , t 


Xk  G Tiki 
Yk  G Ski 


k — 1 ^ ^ t 


Define 
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F0(X0,Y0)  = j 

Then  (P2)  can  be  solved  by  dynamic 
tions: 


0 if  X0  = 0 and  Y0  = 0 
oo  otherwise. 

programming  with  the  following  recursive  equa- 


and 


Ft{XuYt)  = Ht(Xt)  + min  {Ft-\(Xt-\,  Yt-\)  + Pt(xt,  Yt-\,yt) 


Xt  = Xt-1  + xt,  Yt  - g(Yt-i,yt), 


Vt  £ Xt,  Yt-i  £ St- 1,  Xt- 1 £ Ft- 1, 

xt  £ Xt(Yt~i,yt)},  V t 


Z — min{  Ft(Dt,  Yt)  \ Yf  £ St  }. 

Note  that  conventional  dynamic  programming  algorithms  treat  the  state  and 
control  variables,  which  are  Xt,  Yt,  Xt,  and  yt,  as  discrete  variables.  Our  approach 
treats  Xt  as  well  as  xt  as  continuous  variables.  Also,  all  the  cost  functions  are  assumed 
to  be  piecewise  linear  functions  with  a finite  number  of  segments. 

Since  all  given  cost  functions  are  piecewise  linear,  it  is  reasonable  to  assume  that 
Ft(Xt,  Yt)  is  also  a piecewise  linear  function  for  a given  Yt  and  for  any  t.  The  following 
theorem  will  be  proven  later. 

Theorem  2 For  any  t and  Yt,  Ft(Xt,Yt)  is  a piecewise  linear  function  of  Xt  with 
finitely  many  pieces. 

Note  that  the  decision  of  yt  and  Yt  is  independent  of  xt  and  Xt  in  (3.8).  Our 

approach  takes  advantage  of  this  property  and  the  discrete  and  finite  property  of  yt 
and  Yt. 

Given  an  Yt,  let  Tt(Yt ) be  the  set  of  all  possible  combinations  of  (Yt-i,yt)  such 
that  g(Yt-i,yt)  = Yt,  that  is, 


Tt(Yt)  — { (Yt- 1,  yt)  | Yt  — g(Yt-i,  yt),  Yt-\  £ St,  yt  £ }, 
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then  (3.8)  can  be  written  as 

Ft{Xt,Yt)  = Ht(Xt)  + min  {Ft-i(Xt-h  Tt_i)  + Pt(xt,  T)_i,  yt)  \ 

Xt  = Xt~\  + Xt,  Xt- 1 € fct-i, 

xt  e Xt{Yt-i,yt),  (Yt-i,Vt)  € Tt(Yt)}.  (3.9) 

for  any  t — 1, . . . , T.  Now  let 

FtiY,-uyt](Xt,Yt)  = min{i?l_1(X«_a,yt_1)  + Pt(*t,yt_1,yt)  | 

Xt  — Xt~ i + Xt_i  G 

xt  G },  (3.10) 

therefore  (3.9)  can  be  expressed  as 

Ft(Xt,  Yt)  = Ht(Xt)  + min  {F^Yt~uyt) (Xt,Yt)  \ 

(Yt-hyt)  G Tt(Yt)}.  (3-11) 

for  any  t — 1, . . . , T.  Note  that  we  have  assumed  the  function  value  of  F}Yt~1,yt\xt,  Yt) 
to  be  infinite  if  Xt  is  out  of  its  domain.  When  a pair  of  (Yt-i->yt)  is  given, 

Fpt~1,yt)(Xt,Yt)  is  similar  to  Ft(Xt)  in  Section  3.1.  Therefore,  Fpr*~1’y,\xt,Yt)  can 
be  evaluated  in  the  same  way  as  Ft(Xt). 

Lemma  3 FtYt~uVt\Xt,Yt)  is  a piecewise  linear  function  of  Xt  with  finitely  many 
pieces  for  any  given  Yt,  Yt _ i , yt,  and  t. 

Proof:  This  lemma  can  be  proven  by  the  similar  arguments  in  the  proof  of  Theo- 

rem 1 in  Section  3.1.  | 

Proof  of  Theorem  2 The  theorem  will  be  proven  by  mathematical  induction. 
When  t = 0,  F0(X0,  To)  is  a piecewise  linear  function  with  one  segment  (the  original 
point)  by  definition.  Now  assume  that  Ft-\(Xt~\,  if_i)  for  any  T)_i  is  a piecewise 
linear  function  with  finite  segments  for  t > 1.  From  Lemma  3,  it  is  easy  to  see 


50 


that  /'( ( Xt , Yt)  in  (3.11)  is  the  addition  of  Ht(Xt ) and  the  envelope  of  several  (finite) 
piecewise  linear  functions  with  finite  segments.  Since  Ht(Xt ) is  assumed  to  be  a 
piecewise  linear  function  with  finite  segments,  Ft(Xt , Yt)  is  also  a piecewise  linear 
function  of  Xt  with  finite  segments.  This  proves  the  theorem.  | 

The  procedure  to  compute  the  objective  value  function  Ft(Xt,Yt)  in  period  t is 
as  follows.  First,  since  Yt  is  discrete,  we  compute  several  continuous  one-dimensional 
functions  corresponding  to  each  value  of  Yt,  respectively.  For  example,  in  the  problem 
(Ml),  we  compute  Ft(Xt,  0)  and  Ft(Xt,  1).  Now  given  a Yt , we  determine  all  possible 
combinations  of  (Yt-i,yt)  for  Yt,  that  is,  the  set  of  Tt(Yt).  For  instance,  in  the  same 
problem,  these  are 


T,(  0)  = {(0.0), (1,0)} 

T,(  1)  = {(0,1), (1,1)}. 

Then  for  each  {Yt-uy,}  in  Tt(Yt).  we  calculate  the  function  in  (3.10).  The  final  step 
is  to  evaluate  the  lower  envelope  of  these  functions  and  then  add  Ht(Xf)  as  expressed 
in  (3.11). 

3.2.3  Models  with  Startup  and  Reservation  Costs 

The  dynamic  lot  size  model  with  startup  and  reservation  costs  is  first  discussed 
in  Karmarkar  and  Schrage  (1985),  where  it  is  a subproblem  of  the  deterministic  mul- 
tiple product  cycling  problem  in  the  Lagrangean  relaxation  scheme.  In  particular, 
the  startup  cost  is  identical  to  the  cost  of  changeover  between  different  items  in  the 
multiple  product  problem.  The  reservation  cost  is  the  opportunity  cost  which  is  the 
value  of  the  corresponding  Lagrangean  multiplier.  Later,  Karmarkar  et  al.  (1987) 
studied  this  model  in  detail  and  provided  an  0(T 2)  dynamic  programming  algorithm 
for  the  uncapacitated  cases  and  a branch-and-bound  algorithm  for  the  capacitated 
cases.  For  the  latter  cases  their  method  performed  very  well  on  small  problems  (6 
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to  12  periods).  For  these  cases,  Sandbothe  (1992)  provided  a forward  dynamic  pro- 
gramming approach  which  uses  heuristic  rules  to  detect  empirical  decision  horizons. 
1 his  algorithm  can  solve  problems  with  20  periods  in  reasonable  cpu  time.  Although 

it  succeeded  in  finding  optimal  solutions  for  all  28  test  problems,  this  algorithm  does 
not  guarantee  finding  an  optimal  solution. 

In  this  model,  there  are  two  discrete  states:  machine  off  (Yt  = 0),  or  machine 
on  (Yt  = 1).  Thus,  we  have  to  maintain  and  update  two  objective  value  functions 
Ft(Xt, 0)  and  Ft(Xt,  1).  Since  these  objective  value  functions  are  same  as  those  in 
the  ordinary  lot  size  models  (with  single  control  variable),  it  is  reasonable  that  the 
computational  effort  of  this  model  is  twice  as  much  as  those  in  Section  3.1.  Therefore, 
the  computational  complexity  of  our  algorithm  is  0(T  log  T)  for  the  uncapacitated 
cases.  Moreover,  our  algorithm  can  be  applied  for  the  capacitated  cases. 

3.2.4  Models  with  Learning  in  Setups 

Early  works  of  the  dynamic  lot  size  model  with  learning  in  setups  are  in  Karwan 
et  al.  (1988),  Chand  (1989)  and  Chand  and  Sethi  (1990).  All  of  these  papers  address 
the  uncapacitated  case.  The  first  two  papers  also  assumed  the  demand  to  be  constant. 
The  last  one  gives  a dynamic  programming  algorithm  with  computational  complexity 
0(T3)  for  the  problems  with  the  total  setup  cost  being  a concave  function  of  the  total 
number  of  setups.  Malik  and  Wang  (1990)  and  Van  Hoesel  and  Wagelmans  (1991) 
took  advantage  of  the  new  algorithm  for  the  uncapacitated  DLSM  by  Wagelmans  et  al. 
(1992).  The  former  solves  a sequence  of  the  uncapacitated  DLSM  problems  in  which 
the  total  number  of  setups  is  fixe  a certain  number,  respectively,  it  also  assumes  that 
the  total  setup  cost  is  a concave  function  of  the  total  number  of  setups.  The  latter 
solves  a sequence  of  the  uncapacitated  DLSM  problems  with  parameterized  setup 
costs.  Both  assume  speculative  motive  is  not  allowed,  that  is,  the  unit  production 
cost  plus  unit  holding  cost  in  the  current  period  is  always  greater  than  or  equal 
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to  the  unit  production  costs  in  later  periods.  Both  algorithms  have  computational 
complexity  of  0(T2). 

In  these  models,  there  are  at  most  T discrete  states.  Thus,  the  computational 
complexity  of  our  algorithm  is  0(T2logT)  for  the  general  uncapacitated  model  and 
0(T 2)  for  the  uncapacitated  model  without  speculative  motive.  In  addition,  there 
is  no  assumption  about  the  pattern  of  the  setup  costs  with  our  algorithm.  Also  our 
algorithm  is  the  first  to  solve  capacitated  models. 

3.2.5  Models  with  Machine  Replacement 

We  combine  the  dynamic  lot  size  model  and  the  machine  replacement  problem. 
We  call  the  new  model  the  dynamic  lot  size  model  with  machine  replacement.  In  this 
model,  we  consider  inefficiencies  which  are  incurred  by  the  machine’s  age.  These  inef- 
ficiencies are  reflected  by  an  increasing  production  cost  and  a decreasing  production 
capacity. 

The  machine’s  age  can  be  defined  as  the  number  of  setups  as  we  do  in  the 
previous  section.  It  can  also  be  defined  according  to  the  cumulative  production  level. 
For  example,  the  age  increases  by  one  after  producing  a thousand  units  of  product. 
The  only  action  which  concerns  the  machine’s  age  is  to  replace  the  machine  in  our 
preliminary  model.  Actually,  it  is  possible  to  consider  the  action  to  upgrade  the 
machine  to  a certain  age  instead  of  replacing  it  with  a brand  new  machine. 

The  computational  effort  for  this  model  is  small  with  our  algorithm.  In  T-period 
problems,  it  will  take  T times  the  computation  time  of  the  ordinary  lot  size  model. 
When  T — 96  for  example,  it  will  take  only  about  133  seconds  on  average  on  a Vax 
8600  according  to  our  computational  results  in  Section  3.1. 

3.2.6  A Numerical  Example 

The  following  three-period  problem  is  served  as  an  illustration  of  the  algorithm. 
This  is  a dynamic  lot  size  model  with  reservation  and  startup  costs.  In  this  example, 
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rt  is  the  reservation  cost,  and  St  is  startup  cost.  No  holding  cost  is  considered  in  this 
example. 


t 

l 

2 

3 

df 

10 

10 

30 

n 

10 

10 

10 

st 

40 

40 

40 

Pt 

2 

3 

1 

The  optimal  value  functions  are  represented  in  Figure  3.2(a)  to  Figure  3.2(f). 
Figure  3.2(a),  Figure  3.2(c)  and  Figure  3.2(e)  are  the  cases  of  Yt  = 0,  that  is,  machine 
is  off  in  the  corresponding  period.  On  the  other  hand,  Figure  3.2(b),  Figure  3.2(d)  and 
b igure  3.2(f)  are  the  cases  of  Yt  = 1,  that  is,  the  machine  is  on  in  the  corresponding 
period. 

The  machine  can  not  be  off  in  period  1,  thus  the  production  cost  in  Figure  3.2(a) 
is  infinite.  Pp^  in  Figure  3.2(b)  is  the  production  line  of  period  1,  which  is  equal  to 
40  + 10  + 2Xi.  Obviously,  Fi(Xi,  Y\  = 1)  = Pp\ 

P2  ^ in  Figure  3.2(c)  is  a copy  of  Fi(Xi,Yi  — 1),  which  means  no  production  in 
this  period.  P2  in  Figure  3.2(d)  is  a result  of  paying  the  reservation  cost  to  keep  the 

machine  on  in  period  2.  P2(2)  equals  Pa(1)  + 10.  P2(3)  is  the  production  line  in  period  3, 
which  is 

P2(3)  = Fi(10,  Yx  = 1)  + 10  + 3(V2  - 10)  = 80  + 3{X2  - 10). 

F2{X2,Y2  = 1)  is  the  lower  envelope  of  P22^  and  P23\ 

P3(1)  and  P3(2)  in  Figure  3.2(e)  is  a copy  of  F2(X2,Y2  = 0)  and  F2(X2,Y2  = 
1),  respectively.  P3(3)  and  P3(4)  in  Figure  3.2(f)  is  the  production  line  set  up  from 
F2(X2,Y2  — 1)  and  F2(X2,Y2  = 0),  respectively.  These  production  lines  are 

P3(3)  = F2(20,  Y2  = 1)  + 10  + (X3  - 20)  = 110  + (X3  - 20), 

P3(4)  = F2(20,  Y2  = 0)  + 10  + 40  + (X3  - 20)  = 140  + (X3  - 20). 

F3(x3,  Y3  = 0)  is  the  lower  envelope  of  P3(1)  and  P3(2),  and  F3(X3,  Y3  = 1)  is  the  lower 
envelope  of  pj3^  and  P^\ 
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The  optimal  objective  value  of  this  example  is  the  minimal  value  of  F3(X3  = 
50,  Y3  = 0)  and  ^3(^3  = 50, 13  = 1).  The  optimal  value  is: 

Z*  = min{F3pf3  = 50,  Y3  = 0),  F3(X3  = 50,  Y3  = 1)}  = min{140, 150}  - 140. 

The  optimal  plan  is  x\  = 20,  x2  = 0,  x3  = 30  and  the  machine  is  kept  on  in  Period  2. 

3.3  Error  Bound 

The  dynamic  lot  size  model  presumes  that  the  data,  costs  and  demands,  are 
known  precisely.  In  the  real  world,  these  are  obtained  through  forecasting,  and  there- 
fore, reliable  data  are  usually  only  available  for  the  first  several  periods.  This  reliable 
data  horizon  may  not  be  long  enough,  in  that  information  beyond  this  horizon  may 
affect  decisions  in  the  first  or  first  few  periods.  In  the  literature  there  are  two  ap- 
proaches to  deal  with  this  insufficiency  of  information.  The  first  approach  is  to  extend 
the  data  horizon.  The  second  one  is  to  find  a worst-case  bound  on  the  error  induced 
by  imposing  a finite  data  horizon  on  the  model. 

Before  reviewing  the  literature,  we  define  some  terminology  related  to  these  two 
approaches.  Following  Denardo  and  Lee  (1991),  the  number  of  periods,  T,  in  which 
the  needed  data  are  reliable  is  called  a data  horizon  (some  authors  called  it  a time 
horizon ).  On  the  other  hand,  the  number  of  periods  for  which  planning  must  occur, 
TV,  is  called  the  natural  horizon.  Since  we  have  to  determine  the  production  plan  for 
the  first  period  or  the  first  few  periods  based  on  the  data  over  only  the  data  horizon, 
the  question  is:  how  much  error  occurs  if  this  plan  is  implemented?  The  error  bound 
is  the  maximal  possible  value  of  this  error.  The  number  of  these  initial  few  periods 
in  which  the  production  levels  are  determined  and  implemented  constitute  a decision 
horizon  (or  a planning  horizon) . A data  horizon  is  also  called  a forecast  horizon  if  the 
future  information  after  the  data  horizon  does  not  affect  this  initial  production  plan. 
In  other  words,  a data  horizon  is  a forecast  horizon  if  the  error  bound  is  zero.  The 
minimal  forecast  horizon  is  called  the  perfect  horizon. 
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Figure  3.2. 
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Several  authors  have  provided  algorithms  to  find  a forecast  horizon,  if  it  exists. 
We  cite  Wagner  and  Whitin  (1958),  Zabel  (1964),  Eppen  et  al.  (1969),  Blackburn  and 
Kunreuther  (1974),  Lundin  and  Morton  (1975),  Bensoussan  et.  al.  (1983),  Chand 
and  Morton  (1986),  Bean  et.  al.  (1987),  Chand  et.  al.  (1990),  and  Federgruen  and 
Tzur  (1990).  Most  literature  assumes  that  speculative  motive  is  not  allowed,  namely, 
Pt  + ht  > pt+i  for  all  t,  where  pt  and  ht  are  the  unit  production  cost  and  unit  holding 
cost  at  period  t , respectively.  Since  the  unit  production  cost  plus  the  unit  holding  cost 
is  not  less  than  the  unit  production  cost  of  any  later  period,  there  is  no  (speculative) 
motive  to  produce  earlier.  The  only  exceptions  are  works  by  Eppen  et.  al.  (1969), 
Bean  et  al.  (1987),  and  Federgruen  and  Tzur  (1990).  Furthermore,  Bean  et  al.  (1987) 
deals  with  discounted  dynamic  lot  size  models.  Bean  et  al.  (1984)  and  Chand  et  al. 
(1990)  studied  conditions  for  the  existence  of  a forecast  horizon.  Chand  and  Morton 

(1986),  and  Federgruen  and  Tzur  (1990)  proposed  algorithms  to  find  the  minimal 
forecast  horizon,  if  it  exists. 

Lee  and  Denardo  (1986)  and  Denardo  and  Lee  (1991)  proposed  methods  to  find 
the  worst-case  error  bound  for  fixing  the  time  horizon  when  speculative  motive  is  not 
allowed.  Federgruen  and  Tzur  (1990),  and  Chen  and  Lee  (1991)  provided  algorithms 
to  find  an  error  bound  for  the  problem  with  a general  cost  structure.  The  former 
finds  a “upper  bound”,  and  the  latter  finds  a “tight  bound”. 

Unlike  the  works  in  the  literature  which  compute  the  error  bound  after  the 
production  plan  of  the  first  few  periods  is  determined,  the  approach  of  this  paper 
decides  the  plan  for  which  the  corresponding  error  bound  is  the  smallest. 

3.3.1  The  Dynamic  Lot  Size  Model  and  the  Error  Bound 

In  this  section,  we  first  define  the  mathematical  formulation  for  the  single-item 
uncapacitated  dynamic  lot  size  model.  We  also  define  mathematically  the  error  bound 
for  this  dynamic  lot  size  model. 
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Formulation 


Let  N be  the  length  of  the  natural  horizon.  The  problem  data  consists  of  dt , 
Kt,  Pt  and  ht , which  denote  the  demand,  setup  cost,  unit  production  cost  and  unit 
holding  cost  for  period  £,  respectively.  Assume  that  backlogging  is  not  allowed  and 
all  the  costs  are  non-negative.  Then  the  dynamic  lot  size  model  can  be  formulated 
as: 


N 

(PO)  Z*=  min  V (Kt8(xt)  + ptxt  + ht(Xt  - Dt)) 

yXh...1XN)  ^ 2 

subject  to 

t t 

* = £*<>£*  = A,  t = l,...,N 

!=1  i=l 

xt  > 0,  t = 1, . . . , N 


where  xt  the  production  level  in  period  t,  and  8(xt)  = 0 when  xt  = 0,  otherwise 
8(xt)  — 1,  for  any  t. 

Note  that  (P0)  can  be  easily  transformed  into  an  equivalent  problem  without 

holding  costs.  Replacing  Xt  in  the  objective  function  with  xx  H \-xt,  the  objective 

function  of  (P0)  becomes 


N N 

min  ( Kt8(xt ) + ctxt)  — y]  htDt 


where  ct  = pt  + ht  + • • • + for  any  t = 1, . . . , N.  Therefore,  when  the  time  horizon 
is  fixed,  a dynamic  lot  size  problem  can  be  transformed  into  an  equivalent  problem 
without  holding  costs  in  which  ct  are  the  unit  production  costs. 

Let  Ft(Xt)  be  the  optimal  value  function  of  (P0)  considering  only  the  first  t 
periods  with  Xt  as  the  final  cumulative  production  level,  that  is, 

t 

Ft(Xt)  = min  J2  (Kk$(xk)  + PkXk  + hk(Xk  - Dk)) 


subject  to 


Xk  = J2xi > J2di  = Dk,  fc  = i,...,<-i 

t=i  t=i 

Xk  > 0,  k = l,...,t. 
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Then, 

Z*  = Fjv(Z)jv)- 

A recursive  relation  for  Ft(Xt ) is  presented  below.  For  notational  convenience, 
the  holding  costs  are  assumed  to  be  zero.  Define 

p (Y]  _ / 0 for  X0  = 0 

\ oo  for  X0  > 0, 

then  Ft(Xt),  t > 1,  could  be  calculated  by  the  following  recursive  formulation: 

Ft(Xt)  = min{f7i_i(Aii),  Pt{X)}  for  Xt  > Dt 

where 


Pt(Xt)  — Ft-i(Di-i)  + Kt  + ct(Xt  — Dt^ i)  for  Xt  > Dt. 

We  call  Pt(Xt)  the  production  line  of  period  t. 

Chen  and  Lee  (1991)  have  shown  that  Ft(Xt)  is  a piecewise  linear  and  non- 
decreasing concave  function  of  Xt.  It  can  be  obtained  by  a dynamic  programming 
algorithm  with  complexity  O(t\ogt).  Further,  there  are  at  most  t pieces  in  Ft(Xt) 
for  any  t < N. 

A Numerical  Example 

The  following  three-period  example  is  used  to  illustrate  the  above  results.  Note 
that  there  is  no  holding  cost  in  this  example. 


t 
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dt 
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It  is  obvious  that  Pi(Xi)  = 10  + 5AX  (Px  in  Figure  3.3)  and  F^Xi)  = Pi{Xi)  over 
X\  > 10.  When  t = 2,  P2(X2)  is  constructed  based  on  point  a,  that  is,  the  demand 
before  period  2 must  be  satisfied  by  the  production  in  period  1.  It  is  easy  to  see  that 
P2{X2)  = Fi(10)  + 40  + 3(X2  — 10)  = 60  + 40  + 3(^3  — 10)  and  F2(X2)  is  the  lower 
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Figure  3.3. 


envelope  of  P1  and  P2  over  X2  > 20.  Similarly,  P3  = F2( 20)  + 5 + 4(X3  — 20)  = 

110  + 5 + 4(X3  — 20)  and  F3(X3)  is  the  lower  envelope  of  Pi,  P2  and  P3  over  X3  > 30, 
that  is, 


F3(X3)  = 


110  T 5 + 4(X3  — 20)  for  30  < X3  < 35 
60  + 40  + 3(X3  - 10)  for  35  < X3. 


The  total  demand  D3  = 30  is  located  on  P3,  thus  the  last  period  with  a positive 
production  level  is  Period  3.  The  production  of  Period  3 is  based  on  Fi(20).  This 

implies  that  = 20  and  x2  = 0.  Therefore,  the  optimal  solution  is  x\  = 20,  x2  = 0, 
and  x3  = 10. 


The  Error  Bound 


Suppose  that  the  decision  horizon  is  k,  that  is,  the  production  levels  in  the  first 
k periods,  xi,...,xk,  are  determined  and  implemented.  Usually,  k is  smaller  than  the 
data  horizon  T . For  notational  convenience,  let  xk  represent  this  production  plan, 
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that  is,  xk  — (a?i, . . . , Xk).  Then,  we  have  the  following  restricted  version  of  (P0): 


N 


(Pk)  Z*(x)=  min  2 (Kts(xt)  + Ptxt  + ht(Xt  - Dt)) 


subject  to 


Xt  = T,Xi^Y,di=  A,  t = l,...,N 

i= 1 i=l 

^ 0,  t = k -L  1, . . . , AT. 


(Note  that  the  first  k terms  of  the  objective  are  constant.)  Similar  to  the  def- 
initions above,  for  t > k,  let  Ft(x  \xt)  be  the  optimal  value  function  of  (Pk)  where 
Xt  is  the  accumulative  production  level,  and  let 

z*i*k)  = F{/\dn). 


We  are  interested  in  the  difference  between  the  optimal  values  of  (Pk)  and 
(P0),  that  is,  Z*P  1 — Z*.  This  error  could  be  calculated  if  all  data  over  the  natural 
horizon  were  known.  Since  forecasting  all  data  over  the  natural  horizon  is  considered 
unrealistic,  we  are  to  find  the  bound  of  this  error  under  the  known  data,  that  is,  over 
the  data  horizon,  T . Let  £t  be  an  instance  of  possible  data  in  period  t,  that  is, 

6 = {Kt,Pt,ht,dt}. 


Thus,  the  error  bound , given  the  first  k production  levels  xk  and  the  data  horizon  T, 
is  defined  as  follow: 


E(xk,T)  = sup  {Z*^-Z*}. 

Our  analysis  in  the  next  section  uses  the  results  of  Theorems  3 and  4 below. 
Similar  results  have  been  shown  in  Chen  and  Lee  (1991).  A new  proof  is  given  because 
of  the  difference  in  the  definition  of  E(xk,  T ).  In  the  literature  as  well  as  in  Chen  and 

Lee  (1991),  the  choice  of  (xi, . . . , xfc)  is  under  the  restriction  of  xx  + • • • + Xk  = Dk , 
that  is,  Ik  = 0. 
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The  first  result  indicates  that  this  error  bound  is  the  maximal  difference  between 


Ft  {Xt)  and  Ft(Xt)  over  all  Xj  > D?.  In  other  words,  it  is  independent  of  the 
future  information  beyond  the  data  horizon  T . This  result  is  not  a surprise  since  only 
the  outgoing  inventory  in  period  T (or  the  cumulative  production  level  Xt)  affects 
the  production  plan  beyond  period  T . That  is,  given  an  Xt  and  an  instance  of  future 
data  (£x-)-i,  • • • ,£n)i  the  optimal  production  plan  (a?x+i,  • • • , #;v)  depends  only  on  Xt-, 
not  on  production  levels  (xi, . . . ,xj).  Hence  for  any  Xt  > Dt , 


Z*^-Z*  = Fi?k)(XT)  - Ft(Xt). 


Theorem  3 


E(xk,T)  = sup  {F^  \XT)  - Ft(Xt)}. 

X.j'  j' 


Lemma  i For  any  (£t+i,  • • • , £n),  and  t e [k  + l,,N], 

sup  {Ff\xt)  - Ft(Xt)}  < sup  - Ff—\ {Xf—i ) } 

*t>Dt 


Proof:  Given  an  either  X * = Dt- 1 such  that 


Ft(Xt)  — Ft-i(X  ) + A t ct(Xt  — X*), 
Ft{xk\xt)  < Ft*l\x*)  + Kt  + Ct(Xt  — X*). 

or  X*  = Xt  such  that 


F{Xt)  = Ft_i(X*), 

Ff]{Xt ) < Fl*\\x*). 

In  both  case,  it  is  true  that 

Ff\xt)  - Ft(Xt)  < Ft(lk\x*)  - Ft^{X*) 

< sup  {F^ {Xt.t)  - F^iXt.,)} . 
Xt-i>Dt-i 
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Since  the  above  inequality  holds  for  every  X t > Dt  and  the  right  hand  side  is  inde- 
pendent to  Xt,  the  following  inequality  holds: 


sup  {rffc)(xf)  - Ft(xt)}  < 

Xt>Dt 


sup  {Ft{xk)(Xt- 1)  - Ft(Xt-i)}. 


I 

Proof  of  Theorem  3:  We  first  show  that  E(xk,T ) > sup {F^k\XT)  - FT(XT)  : 

Xt  > Dt}.  Let  X*  > Dt  be  such  that  Fj  *(V*)  — FT(X*)  = sup{FyX*'(XT)  — 
Ft(Xt)  : Xt  > Dt}.  Let  DN  equal  to  X*,  and  let  I<t  for  all  t > T be  very  large 
numbers  such  that  Ff$k\x*)  = F^\(X*)  = ■■■  = F^k\X*).  Then 


Since 


we  have 


Z*  = Fn(Dn)  < Ft{Dn ) = Ft(X*), 


= ( Xk ) 


E(x*,T)  = sup  {Z 

> f(/\dn)-fn(dn) 


> Fi?k)(X*)-FT{X*) 

— sup  {Fj.  ^(Vr)  — F t ( X'i ) } 

Xt  >Dt 


z*} 


(3.12) 


Now,  we  show  that  E(xk,T)  < supfi^^Xr)  - FT(XT)  : XT  > DT}.  From 
Lemma  4,  it  is  straightforward  that  for  any  (£t+ i,  • • • ,6v), 


Z 


(**) 


N 


Zn  < 


< 


sup  {Fj^k\xN)  - FN(XN)} 

xn>dn 

sup  {F^\(Xff-\)  - Fjv-i(Vjv-i)} 


< sup  {F^x  \Xt)  - Ft(Xt)}. 

Xt>  Dt 
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Note  that  the  above  inequality  holds  for  any  (£r+i,...,£/ v)  and  the  right  side  is 
independent  to  any  future  data.  Therefore 


E(xk,T)  = 


< 


sup  {Z(^k)  - ZN } 

»£at) 

sup  {Ft  V-Xx)  — Ft(Xt){ 
xt>dt 


(3.13) 


From  (3.12)  and  (3.13),  we  have 


E(x\T)  = sup  {F^xk) (XT)  - Ft(Xt)} . 

xt>dt 


Since  Ft  (Xt)  and  Ft(Xt)  are  piecewise  linear  concave  functions,  it  is  easy  to 
see  that  the  maximum  of  FT  {Xt)  — Ft{Xt)  occurs  at  one  of  the  break  points  in 

Ff](X  T). 


Theorem  4 Let  s be  the  period  such  that  cs  = ps  + hs  H 1 - hT  is  the  smallest 

slope  in  Ft(Xt),  then 


E{xk,T)  = 


_ Ft(Xt)}  ifs>k 

OO  otherwise, 


where  B*{T)  is  the  set  of  Xt  at  which  Xt  Z Dt  and  Fj'  \Xt)  changes  its  slope. 


Note  that  there  are  at  most  T break  points  in  Fjf  Thus,  for  the  data 

horizon  T,  given  an  instance  xk,  the  error  bound  E(xk,  T ) can  be  evaluated  with  effort 
0(7  log  T)  and  0(T ) for  the  cases  without  and  with  speculative  motive,  respectively, 


Numerical  Example  (continued) 

Suppose  that  the  decision  horizon  is  one.  The  conventional  approach  would  be 
to  set  the  production  level  of  the  first  period  according  to  the  optimal  solution  of  the 
three-period  problem,  that  is,  x1  = 20,  and  to  compute  F^Xl=20)(X 3).  If  this  is  done, 
point  a (aq  = 10)  in  Figure  3.3  is  not  available.  Actually,  only  point  b on  Pi(Na) 
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is  available  in  this  restricted  case.  Therefore,  P2(X2)  is  replaced  by  P'PG)  in  the 
restricted  model,  that  is,  P£{X2)  = Fi(20)  + 40  + 3(X2  - 20)  = 110  + 40  + 3(X3  - 20) 
(see  Figure  3.4).  F^Xi~20\x3)  is  the  lower  envelope  of  Pi,  P'2  and  P3  over  X3  > 30: 

p(Xl=20)(x  , = f 110  + 5 + 4{X3  - 20)  for  30  < X3  < 55 
3 1 3>  \ 110  + 40  + 3(^3-20)  for  55  < X3. 

F3(X3)  and  Ft=20)(X3)  are  both  in  Figure  3.5.  Note  that  these  two  functions  are 
parallel  to  each  other  when  X3  > 55.  The  error  bound,  which  occurs  at  point  c 
(X3  = 55),  is 

E{Xl  = 20,  T = 3)  = Ff1=20)(55)  - F3(55)  = 20. 


3.3.2  Minimizing  the  Error  Bound 

As  defined  above,  the  error  bound  E(xk,  T)  is  a function  of  xk  and  T.  Thus, 
this  error  bound  can  be  reduced  by  varying  xk  or  T.  It  has  been  mentioned  in  Lee 
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Figure  3.5. 

and  Denardo  (1986)  that  the  error  bound  can  be  reduced  by  increasing  T.  In  this 
section,  we  discuss  the  method  to  minimize  the  error  bound  by  selecting  an  . 

To  simplify  this  problem,  we  first  discuss  the  case  of  Ar  = 1 . In  other  words, 
given  a data  horizon  of  T periods,  we  are  to  find  the  minimal  error  bound  over  all 
possible  X\\ 

m 

E'(T)  = inf  {£(*,,  T)}. 

The  case  of  k > 1 will  be  discussed  later. 

Let  St  = [Dt,  Dt+i),  for  t = 1, . . . , T — 1,  and  St  = [Dt,  oo). 

Theorem  5 For  fixed  T , if  Pl  + hx  + • • • + > Pi  for  all  2 < i < T,  then  E(x i , T) 

is  an  increasing  function  of  xx  for  xx  e St,  t = 1, . . . , T. 

Proof:  Let  a,  6 6 St  and  a < b,  we  first  show  that  Fp=b)(XT)  > F^Xl=a)(XT), 

for  Xj  > Dt . Let  (b,  0, ... , 0,  = b\  xt-+1, . . . , xt)  be  the  corresponding  production 
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plan  for  Tj  1 \Xt)  where  2 < i < T.  Let  Z be  the  objective  value  of  the  plan 
(a,  0, . . . , 0,  X{  = b +6  — a,  . . . , xt)>  Then,  1 ^ {<Xrp')  > Z since  Pi  + h\  + • • • + 

hi- 1 ^ P*-  Also,  it  is  easy  to  see  that  Z > Therefore  > 

Z > F^1=a)(XT). 

Let  X£  > be  such  that 

4',=,)(Xf)  - FrW)  = sup  {F<*'=“>(XT)  - Ft(A:t)}, 

Xj'^>D  j' 

then 

E(b,T)  - E(a,T)  = sup  {F^1~b\xT)  — FT(XT)} 

xt>dt 

- sup  {Fj1~<1\Xt)  - Ft(Xt)} 
xt>dt 

> F^Xl=b\x*T)  - Ft{X*t)  - FjXl=a\x^)  + Ft(Xj) 

> 0. 

Hence,  E(xi,T)  is  an  increasing  function  of  xx  for  xx  € St}  t = 1, . . . , T.  | 

Corollary  1 Given  a T , if  px  + hx  + • • • -f-  /i;_ x > px  for  all  2 < i < T , then 

E*(T)=  min  {E(xuT)}, 

x\€{Di,—,Dt} 

and  E*(T)  can  be  obtained  with  computational  complexity  of  ()( T2  log  / ) for  the 
general  case  and  0(T2)  when  speculative  motive  is  allowed. 

Proof:  According  to  Theorem  5 and  the  definition  of  E*(T ),  it  is  easy  to  see  that 

E*(T)  = min  {E(x1,T)}. 

Since  given  an  xu  E(xuT)  can  be  obtained  in  O(TlogT)  for  the  general  case 
and  in  0(T)  when  speculative  motive  is  allowed,  by  evaluating  all  possible  xx  6 

. • ,-Dr},  E*(T)  can  be  obtained  in  0(T2\ogT)  or  0(T 2)  for  the  cases  with  or 
without  speculative  motive,  respectively.  | 
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Numerical  Example  (continued) 

The  error  bound  E(x\  =20 , T = 3)  is  20  as  shown  in  the  previous  section.  When 
x\  is  restricted  to  be  10, 

FjXl=10)(X3)  = 60  + 40  + 3(X3  - 10)  for  30  < X3. 
and  the  error  bound,  which  occurs  at  X3  = 30,  is 

E(Xl  = 10,T  = 3)  = F3(xi=10)( 30)  - Fs(30)  = 5. 

Also  when  x-i  is  restricted  to  be  30, 

ri(^i=3o)  / -w-  \ _ f 160  + 5 + 4(X3  — 30)  for  30  < X3  < 65 

3 1 3)  ~ \ 160  + 40  + 3(X3  - 30)  for  65  < X3. 

and  the  error  bound,  which  occurs  at  X3  = 65,  is 

E{Xl  = 30,  T = 3)  = Tf1=30)(65)  - X3(65)  = 40. 

Therefore,  the  minimal  error  bound  E*(T)  = min{5,20,40}  = 5. 

Remark  If  we  fix  X\  = 20  which  is  obtained  by  solving  the  three-period  problem, 
the  error  bound  is  20.  However,  using  xi  = 10  gives  a smaller  error  bound,  which,  by 
Corollary  1,  is  minimal. 

Case  of  k > 1 

In  the  case  of  k > 1,  only  the  summation  of  xk,  that  is,  xa  + j-  Xk,  is  critical. 

Let  S{ Xk)  = {(®i, . . . , Xk)  | #i  + • • • + Xk  = Xk}  and  (®i, . . . , Xk)*  £ S(Xk)  be  one  of 
the  optimal  solutions  for  Fk(Xk).  Then,  it  is  not  difficult  to  prove  that 

E((Xl, . . . ,Xk)*,T)  < E((xi, . . . , xk),  T)  for  any  (®i, . . . ,xk)  £ S(Xk). 

In  other  words,  the  error  bound  can  be  treated  as  a function  of  Xk.  Therefore  given 
a data  horizon  T,  we  are  to  find  the  minimal  error  bound  over  all  possible  Xk: 

E*(T)  = xMD{E(Xk,T)}. 

The  results  in  this  given  are  still  valid  for  this  case. 


CHAPTER  4 

A HEURISTIC  FOR  MULTI-ITEM  DYNAMIC  LOT  SIZE  MODELS 


In  this  chapter,  we  propose  a heuristic  algorithm  for  the  MLSM  which  is  based 
on  the  algorithm  for  single-item  capacitated  lot  size  models.  Experimental  results 
show  that  this  heuristic  algorithm  is  very  effective  in  finding  a good  primal  feasible 
solution  of  the  multi-item  model. 

4.1  Formulation  and  Lagrangian  Relaxation 

We  need  the  following  notation  to  describe  the  general  multi-item  dynamic  lot 
size  model. 

da  = demand  of  item  i in  period  t. 

xn  = production  level  of  item  i in  period  t. 

Dit  = cumulative  demand,  that  is,  Du  = Y?k=i 

Xu  = cumulative  production,  that  is,  Xu  = YX= i xik- 

Xu  = a set  of  feasible  production  levels,  that  is,  xu  (E  Xu- 

Hit  — a set  of  feasible  cumulative  production  levels,  that  is,  Xlt  £ 7 Zu- 

Pit(xit)  — cost  function  for  producing  xxi  units  of  item  i in  period  t. 

Hit(Xu)  = cost  function  for  carrying  Xu  — Du  units  inventory 
of  item  i at  the  end  of  period  t. 

Gu(xu)  — resource  consumed  by  production  of  item  i in  period  t. 
ct  = common  resource  limit  in  period  t. 

Without  loss  of  generality,  we  may  assume  that  Wo  = 0 and  X,T  = At  for  all 
i.  And  for  notational  convenience,  let  Di0  = 0 for  all  i.  The  general  MLSM  can  be 
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formulated  as  follows: 


such  that 


(P3)  Z*  = min^^  (Ptt(xlt)  + Htt{Xlt)) 


* = 1 2—1 

X Xij—\  T X{f^  V i,t 

(4.1) 

i 

Git(xit)<cu  Vi 

(4.2) 

2 — 1 

xit  € V i,2 

(4.3) 

xitenit,  Vi,  t. 

(4.4) 

The  constraints  (4.1),  (4.3),  and  (4.4)  in  the  above  model  are  similar  to  those  in 
the  generic  single-item  dynamic  lot  size  model.  The  only  constraints  in  which  items 
are  related  to  each  other  are  the  common  resource  constraints  (4.2). 

Usually,  constraints  (4.2)  are  relaxed  in  the  multi-item  dynamic  lot  size  model. 
Let  7 rt  be  the  corresponding  multipliers.  The  Lagrangian  dual  is  as  follows: 

T i t 

(D)  L(  tt)  = min  ^ ^ ^ (x(ct) 

t=l  2 — 1 t=l 

such  that 


X it  — X{j—\ 

+ xiti  V i,  t 

(4.5) 

xit  £ Xit, 

Vi,  t 

(4.6) 

Xu  G 7 Zn, 

V i,  t. 

(4.7) 

Then  evaluation  of  L(tt)  can  be  separated  into  / single-item  problems  with  a particular 
7 r given.  If  Pit(xit ),  Hit(Xit),  and  Git(xit ) are  piecewise  linear  functions,  these  single- 
item uncapacitated  problems  can  be  solved  with  any  of  several  algorithms  including 
those  presented  in  the  last  chapter.  Also,  it  is  well  known  that  T(7r)  is  a concave 
function  and  the  dual  problem  is  to  maximize  L(k)  with  7 r non-negative. 
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For  simplification,  we  first  consider  the  problems  without  setup  time.  That  is, 


= KitS(xit)  + pitxit 


Hu  (Xn)  = 

-Dit) 

Git(xn)  — 

xit  = 

{ %it 

0 < xit  } 

Hit  = 

{Xit 

Ditt- 1 < Xn  < Dix  } 

where  Ka , pa,  ha  and  fin  are  the  setup  cost,  unit  production  cost,  unit  holding  cost 
and  resource  consumption  rate  of  item  i in  period  f,  (3it  is  0 if  xit  = 0,  or  1 if  xit  > 0. 

At  the  end  of  this  chapter,  we  discuss  an  algorithm  for  the  model  with  setup 
time,  that  is, 

^it(xit)  ^it^(xit)  + ftitXit 


where  an  is  the  resource  required  if  there  is  set  up  for  item  i in  period  t. 


4.2  Algorithm 

Two  issues  are  of  concern  in  solving  ( D ) as  well  as  (PS).  First,  how  are  the 
multipliers  7Tf  updated  to  obtain  dual  optimality?  Second,  how  can  a (primal)  feasible 
solution  be  recovered  from  a dual  feasible  solution?  Since  L(n)  is  a concave  function 
of  7 r,  there  are  many  techniques  to  maximize  L(ir):  such  as,  subgradient  methods 
(Held  et  al.,  1974),  the  dual  cutting  plane  method,  generalized  linear  programming, 
and  bundle  methods  (Kiwiel,  1983).  Among  the  methods,  a subgradient  method  is 
the  most  frequently  used  method  because  that  it  performs  better. 

In  our  algorithm,  we  will  use  the  subgradient  technique  used  by  Thizy  and  Van 
Wassenhove  (1985).  The  formula  for  updating  multipliers  is 


Trt+1  = max  {0, 7r(  + 74 


T*  - Ujrt) 

lk(7r0ll2 


9(*t)}, 
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where  g(irt)  is  the  subgradient  at  irt  and  L*  is  the  best  primal  feasible  solution  found 
so  far.  The  parameter  7 is  set  to  2 at  the  beginning,  then  it  is  divided  by  2 for  every 
T consecutive  non-improving  iterations. 

In  Thizy  and  Van  Wassenhove  (1985),  after  each  subgradient  iteration,  a trans- 
portation problem  is  solved  in  order  to  obtain  a feasible  solution.  Let  {xit}  be  the 
solution  corresponding  to  a dual  value  7 r.  Then  the  production  period  is  determined 
for  each  period  and  item,  and  according  to  Bowman  (1956),  (P3)  can  be  represented 
as  a transportation  problem  with  |T'|  sources  and  |T|  x |/|  demands.  Let  0nT  be  the 
amount  of  resource  in  period  t used  to  produce  item  i for  the  demand  of  period  r. 
The  transportation  problem  is  as  follows. 


ITT 

min  htrOitT 

i= 1 *=1  T= 1 

such  that 

T 

^OitT  > Pidir,  V Z,T 

t- 1 
I T 

Out  < <h,  v t 

1=1  T = 1 

Out  = 0 if  Xit  = 0 or  t > r,  V i,  t,  r 
Out  > 0,  V i,t,r 


where 


& itr  — 


0 lip  it  + Y^j=t  hi])  if  t < T 

otherwise 


00 


V i,  t , r. 


The  decision  variables  of  this  formulation  and  that  of  (P3)  have  the  following 
relation;  that  is,  0ltT  = 0 if  xit  = 0,  and  xit  - j-  Y^=t  0 itr  for  all  *,  t. 

Note  that  if  {xu}  has  the  same  setup  periods  for  every  item  as  the  optimal 
solution  of  (P 3),  then  the  above  transportation  problem  solution  will  be  an  optimal 
solution  of  (-P3).  However,  given  an  {xa}  the  above  transportation  problem  may  be 
infeasible,  thus  there  is  no  corresponding  feasible  solution  of  (P3).  Therefore,  Thizy 
and  Van  Wassenhove  (1985)  report  the  best  feasible  solution  found  over  all  iterations. 
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Our  approach  also  uses  dual  iterations  to  determine  {xa}.  However  we  develop 
a new  heuristic  algorithm  to  find  a feasible  solution  instead  of  solving  the  above 
transportation  problem.  This  heuristic  is  discussed  in  the  next  section. 

4.3  Heuristic  Method 

There  are  two  stages  in  our  heuristic  algorithm:  (1)  finding  a feasible  solution, 
and  (2)  improving  the  objective  value. 

4.3.1  Finding  a Feasible  Solution 

Given  a solution  which  satisfies  all  the  constraints  of  (P3)  except  (4.2).  Then, 
a period  is  said  to  be  ‘‘exceeded”  if  the  resource  required  exceeds  the  resource  limit. 
An  item  in  an  exceeded  period  is  said  to  be  a “candidate”  if  the  production  level  of 
this  item  is  positive  in  this  period. 

The  major  aim  is  to  disperse  the  resource  requirement  in  an  exceeded  period  to 
other  periods  by  discouraging  production  in  this  period  and  by  encouraging  set  ups 
or  production  in  other  periods.  The  procedure  is  as  follows. 

1.  Pick  a candidate  item  in  the  first  exceeded  period. 

2.  Change  the  capacity  of  this  item  in  each  period. 

3.  Increase  the  unit  production  cost  of  this  item  in  the  first  exceeded  period. 

4.  Decrease  the  setup  costs  of  this  item  in  all  periods  except  the  first  exceeded 
period. 

5.  Solve  a single-item  capacitated  lot  size  model  of  this  candidate  item. 

6.  Update  the  production  level  of  this  candidate  item  and  go  to  Step  1 until  no 
exceeded  periods  exist. 

In  Step  2 the  capacities  are  changed  in  such  a way  that  the  periods  before  the 
first  exceeded  period  remain  non-exceeded  after  the  updating  in  Step  6.  Therefore, 
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resource  is  reserved  for  other  items  in  the  periods  before  the  first  exceeded  period. 
Only  the  remaining  resource  and  the  resource  which  is  already  used  by  the  candidate 
item  is  available  for  the  candidate  item.  On  the  other  hand,  all  the  resource  is  available 
for  the  candidate  item  in  the  first  exceeded  period  and  the  following  periods. 

Ihe  unit  production  costs  and  setup  costs  are  changed  gradually  in  Step  3 and 
Step  4.  At  the  beginning  of  each  iteration,  a changing  rate  is  assigned  for  the  unit 
production  cost  and  the  setup  cost,  respectively.  These  rates  are  increased  after  all 
the  candidate  items  are  picked.  In  this  way,  the  exceeded  resource  requirement  in  the 
first  exceeded  period  can  be  dispersed  to  other  periods. 

The  heuristic  algorithm  for  finding  a feasible  solution  is  summarized  as  follows. 
In  this  algorithm,  {xn}  is  the  current  solution,  St  is  the  remaining  resource  in  period  t 
(St  < 0 means  that  Period  t is  an  exceeded  period),  and  Rp  and  Rs  are  the  basic 
changing  rates  of  the  unit  production  cost  and  the  setup  cost,  respectively. 


Algorithm  for  Finding  a Feasible  Solution 

1.  Let  {ah*}  be  an  initial  solution  which  satisfies  all  the  constraints  in  (P3)  except 
(4.2). 

2.  Let  Sk  = ck  — Yj\=\(PikXik)  for  all  periods  k,  t = 1, . . . , T. 

3.  Let  t be  the  first  exceeded  period;  if  there  is  no  such  period,  stop  (a  feasible 
solution  is  found). 

4.  Let  j = 0 and  repeat  until  St  > 0. 

(a)  Pick  an  item  i for  which  xit  > 0.  If  i is  the  first  item  selected  for  this  t,  let 

j = j + 1. 

(f*)  Let  Ck  Xtk  T Skj j3ik  for  k = 1, . . . , t — 1;  and  let  Ck  = ckj tk  for 
k = t, . . . , T . 

(c)  Let  Pk  = Pik  X (1  + j x Rp)  if  k = t\  otherwise,  let  P 'k  = pik. 

(d)  Let  Kk  = I<ik  if  A:  = t;  otherwise,  let  Kk  = I<lk  x (1  -j  x Rs  x Qik),  where 

Q ik  (Sk  "f*  (3 ikX tk) / ’ Ck . 

(e)  Solve  the  capacitated  lot  size  model  for  item  i using  Ck , Kk,  Pk,  Hik  as  the 
capacity,  setup  cost,  unit  production  cost,  and  holding  cost  for  period  A, 
k = 1, . . . , T.  Let  Xi  be  the  solution. 

(f)  Let  Sk  = Sk  + /3ik(xik  - xik)  and  xik  = xik,  for  k = 1, . . . , T. 

5.  If  t < T,  go  to  Step  2. 
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4.3.2  Improving  the  Objective  Value 

In  this  stage,  the  given  solution  is  feasible,  so  the  goal  is  to  improve  the  objective 
value  by  assigning  unused  resources  in  each  period  to  some  particular  item.  The 
algorithm  is  similar  to  the  one  in  the  last  section  except  there  is  no  changing  of  unit 
production  cost  or  setup  cost. 


Algorithm  for  Improving  the  Objective  Value 

1.  Let  {x{t}  be  the  initial  feasible  solution. 

2.  Let  St  = Ct  — J2i=i(flitXit)  for  all  t. 

3.  Let  k = 0,  and  repeat  until  k = T. 

(a)  Pick  an  item  i in  turn. 

(b)  Let  Ct  = xit  + Sk//3it,  for  t = 1, . . . , T. 

(c)  Solve  the  capacitated  lot  size  model  for  item  i using  Ct,  Ka , pu,  ha  as  the 
capacity,  setup  cost,  unit  production  cost,  and  holding  cost  for  period  £, 
t = 1, . . . , T.  Let  X{  be  the  solution. 

(d)  If  Xi  = Xi,  let  k = k + 1;  else,  let  k — 0. 

(e)  Let  Sk  = Sk  + /3ik{xik  - xik)  and  xik  = xik , for  k = 1, . . . , T. 


4.4  Experimental  Results 

We  compare  the  above  algorithm  with  those  in  Thizy  and  Van  Wassenhove 
(1985)  and  Brannlund  (1992). 

In  Brannlund  (1992),  the  approach  is  similar  to  that  in  Thizy  and  Van  Wassen- 
hove (1985).  The  former  uses  a bundle  technique  to  update  multipliers  instead  of  the 
subgradient  method,  and  the  transportation  problem  is  slightly  modified. 

In  our  proposed  heuristic  algorithm,  the  parameters  Rp  and  Rs  are  set  to  0.05 
and  0.1,  respectively. 

Twelve  test  problems  (T  = 8 and  I = 8)  from  Brannlund  (1992)  are  tested.  The 
results  are  listed  in  Table  4.1.  The  optimal  solutions  (in  Column  2)  and  the  results 
in  Column  3 and  4 are  from  Brannlund  (1992).  Brannlund  only  provided  the  best 
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Table  4.1 


Capacity 

Optimal 

primal 

T&VW 

Brannlund 

Proposed  Method 

Best 
primal 
(50  eval) 

Best 
primal 
(50  eval) 

Primal 
solution 
(1  eval) 

Best 
primal 
(50  eval) 

Best 

found 

at 

(Very  high  setup  costs) 

Tight 

29740 

30490 

30090 

31170 

30380 

12 

Medium- Tight 

26030 

27430 

27060 

26290 

26290 

1 

Medium- Loose 

20920 

20920 

20920 

21020 

21020 

1 

Loose 

19210 

19800 

19800 

19210 

19210 

1 

(High  setup  costs) 

Tight 

8430 

8530 

8530 

8480 

8480 

1 

Medium-Tight 

7910 

8250 

7970 

7970 

7910 

2 

Medium- Loose 

7610 

7610 

7610 

7650 

7610 

31 

Loose 

7520 

7520 

7520 

7520 

7520 

1 

(Low  setup  costs) 

Tight 

2900 

2900 

2900 

2910 

2900 

2 

Medium-Tight 

2900 

2930 

2930 

2925 

2925 

1 

Medium-Loose 

2865 

2865(1) 

2865(1) 

2865 

2865 

1 

Loose 

2865 

2865(1) 

2865(1) 

2865 

2865 

1 

T&VW:  subgradient  technique  + transportation  problem 

(Thizy  and  Van  Wassenhove,  1985) 

Brannlund:  bundle  technique  + transportation  problem  (Brannlund,  1992) 

Proposed  Method:  subgradient  technique  + proposed  heuristic  algorithm 


feasible  solution  in  50  iterations,  and  no  information  about  in  which  iteration  these 


solutions  are  obtained. 

These  results  show  that  our  heuristic  algorithm  is  very  effective.  Even  when 
only  the  pure  heuristic  algorithm  is  used  (see  Column  5),  the  largest  error  is  within 

5%. 


CHAPTER  5 
CONCLUSIONS 


We  have  developed  a dynamic  programming  algorithm  for  general  single-item 
dynamic  lot  size  models  (DLSM).  This  algorithm  performs  better  than  any  existing 
algorithm  for  the  NP-hard  versions  of  the  DLSM.  This  study  also  shows  that  the  new 

algorithm  for  DLSM  can  serve  as  a basic  building  block  for  the  multi-item  dynamic 
lot  size  model  (MLSM). 

Unlike  conventional  dynamic  programming  algorithms,  the  above  algorithm  uti- 
lizes a continuous  state  variable.  By  assuming  the  cost  functions  to  be  piecewise 
linear,  we  have  shown  that  the  optimal  value  functions  are  also  piecewise  linear,  and 
these  optimal  value  functions  can  be  updated  by  evaluating  the  lower  envelope  of 
finitely  many  line  segments.  Extensive  computational  experience  shows  that  in  prac- 
tice the  cpu  time  of  this  algorithm  increases  only  at  an  approximately  cubic  rate  in 
the  number  of  time  periods  of  the  NP-hard  version  of  the  problems. 

This  new  algorithm  can  apply  to  most  variations  of  the  single-item  dynamic 
lot  size  model  in  the  literature.  For  example  as  discussed  in  previous  chapters,  the 
Wagner-Whitin  model,  the  uncapacitated  DLSM,  the  DLSM  with  capacity  on  pro- 
duction,  the  DLSM  with  capacity  on  inventory,  the  DLSM  with  backlogging,  the 
DLSM  with  learning  in  setup,  the  DLSM  with  startup  and  reservation  costs,  and  the 
DLSM  with  machine  replacement  can  all  be  solved  with  this  method. 

We  also  investigate  the  error  bound  issue  of  the  uncapacitated  DLSM.  Given 
a production  plan  (up  to  planning  horizon),  we  provide  an  algorithm,  which  can 
evaluate  the  error  bound  in  0(T\ogT)  where  T is  the  data  horizon.  In  addition,  we 
propose  that  it  is  more  realistic  to  find  a production  plan  such  that  the  possible  error 
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bound  is  minimal  than  to  extend  the  data  horizon  to  lessen  the  error  bound.  We  show 
that  the  minimal  error  bound  can  be  obtained  in  0(T2  log  T)  for  the  uncapacitated 
DLSM.  How  we  can  extend  this  idea  for  the  general  DLSM  is  an  worthy  avenue  for 
future  research. 

We  only  have  a preliminary  study  of  the  MLSM.  However,  the  preliminary  re- 
sults are  very  promising.  Based  on  the  algorithm  for  the  single-item  model,  we  have 
proposed  a heuristic  for  the  MLSM.  Computational  experience  shows  that  it  is  very 
competitive  with  algorithms  in  the  existing  literature  for  the  MLSM  with  a com- 
mon resource.  More  importantly,  this  heuristic  algorithm  can  apply  to  very  general 
MLSMs,  for  instance,  the  MLSM  with  setup  time. 
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